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ABSTRACT 
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the first eigenvector. This dissertation gives necessary and sufficient conditions for the 
case of the real spectrum when the first eigenvector is posittve, and when the first 
eigenvector is nonnegative it gives conditions based on the reducibility of the matrix. For 
the complex spectrum, in the presence ofa positive first eigenvector the problem is shown 


to be equivalent to the irreducible nonnegative inverse eigenvalue problem. 
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I. INTRODUCTION 


This dissertation addresses what appears to be a simple problem: “Given a real, 
square matrix A, when does there exist a nonnegative matrix B such that B™'AB has its 
off-diagonal elements nonnegative?” Based on different bodies of literature, and on the 
application, we could similarly ask when a square linear operator is quasimonotone 
nondecreasing with respect to a cone in the nonnegative orthant, or when a real matrix is 
essentially nonnegative under a nonnegative change of basis. 

This question has its roots im the Perron-Frobenius theory of the early 1900s, and 
as we develop our solution, we trace its history, as well as the history and current state of 
some of its applications. We also discuss similar and related problems to demonstrate how 
this problem fits into the scheme of applied mathematics. 

We present a brief background on cones in Euclidean space and on the theory of 
nonnegative matrices in Chapter II. We use the following convention to present results. 
Theorems (and proofs) from the orginal source have the author and date after the 
theorem number. Ifthe proof comes from another source, we indicate that prior to 
presenting the theorem. For well-known results, we either indicate the source or provide 
our own proof. We indicate the completion of a theorem and proof (or just a theorem for 
results presented without proof) with the (© symbol. 

Chapter HI motivates this question as an unsolved problem from the theory of 
cone-valued vector Lyapunov functions, which can determine stability in dynamical 


systems. We also present some of our results concerning this application. This question 


was first posed for the application of vector Lyapunov functions in 1974, and partial 
answers were given for this setting in 1980 and 1995. 

Chapters IV and V present our onginal solutions to this problem. The problem 
has significantly different solutions depending on whether the spectrum of A is real or 
general, and Chapter IV presents our solution for the real spectrum. In Section IV.A we 
show the necessary condition that the first eigenvector x, of A (associated with the 
greatest eigenvalue A,) be nonnegative, and the sufficient condition that x) be positive. 

The case where x; is nonnegative but not strictly positive has different solutions 
based on the reducibility of the matrix A. Section IV.B addresses the reduced case, where 
we show sufficient conditions when A, is a simple eigenvalue of A, and we have a 
sufficient (but not necessary) condition when A, is not simple. 

In the case where x, is nonnegative and A is irreducible, Section IV.C shows 
that we can use x; to deflate the matrix A and reduce the problem to one of dimension 
n—1. Here, we also give necessary and sufficient conditions for a cone to exist. 

Chapter V presents our solution for the complex spectrum. We show that in the 
presence of a positive first eigenvector, the problem can be reduced to the nonnegative 
inverse eigenvalue problem, a classic unsolved problem from theoretical lmear algebra. 
We present the solution only in terms of irreducible quasimonotone matrices, and we show 
why it is sufficient to consider such matrices. 

Chapter VI presents other direct and related applications from control theory and 
dynamics, as well as from other fields of applied mathematics, and we mention similar 


problems which have been solved or discussed previously from these and other fields. 


Chapter VII summarizes our results, to include the unsolved problems mentioned 
above, and we suggest further research. 

One attractive aspect of our solutions is they are constructive, and hence fairly 
sumple, although there is some numerical instability inherent in some constructions. 
Particularly for applications where the existence of a matrix B is sufficient information for 
the analyst, we present collections of useful necessary and sufficient conditions. We begin 


with some background. 
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Ul. CONES AND NONNEGATIVE MATRICES 


This chapter provides background for the main results. Specifically, we discuss 
cones in Euclidean space and present results from the theory of nonnegative matrices. We 


begin with a discussion of cones. 


A. CONES IN EUCLIDEAN SPACE. 

The employment of cones as subsets of R" began in the 1930s, where along with 
convex polyhedra, they were soon found to be useful in the field of linear optimization and 
in the study of lear mequalities (see, for example, Weyl, 1935; and Goldman and Tucker, 
1956). Dual cones were first employed in 1941 by J. Dieudonné (see Lay, 1982), 
although dual sets to convex sets had been used much earlier (see Helly, 1923). 

The analysis of cones, convex polyhedra, and convex sets in general has become a 
major field of study in its own right, which has its roots in Caratheodory (1907), and 
which can be found, for example, in Sandgren (1954), Karlin (1968), Rockafeller (1970), 
Stoer and Witzgall (1970), Berman (1972), and Lay (1982). We need only the most basic 
definitions and results from those fields, and we begin with some definitions. 

If R" is Euclidean n-space with norm ||-|| and inner product en a subset 
S CR" is convex if and only if x, ye S implies ax+(1-—a)y € S forall a € [0,1]. A 
set K isacone if and onlyif AK CK for A420. A convex cone satisfies both 
definitions, and we require that K= K, where K _ isthe closure of K, so that our cones 
are closed. (It is a result of Farkas (1901) that a set is a closed, convex cone if and only if 


K =(K’),, where K’ is the dual of K, to be defined shortly.) A cone is pointed if and 


only if K 7 (—K) = {0}, and solid if and only if K° is nonempty, where K° is the interior 
of K. An equivalent notion (in Euclidean space) to solid is reproducing, where K is 
reproducing if and only if K- K=R". (it is a result of Krein and Rutman (1948) that a 
cone is pointed if and only if its dual is solid.) 

A convex cone K 1s polyhedral if and only if it is generated by finitely many 
vectors, and if the number of independent vectors is equal to n, then K is simplicial. A 
closed, pointed convex cone is called proper. The result that a proper cone is generated 
by its extremals is a special case of the Krein-Milman Theorem. The only cones we 
consider for our applications are proper, simplicial cones. 

Since a proper, simplicial cone is generated by n independent extremal vectors bi, 


we consider the nonsingular matrix B = [b,,..., bh] € R°", and we denote the cone 


n 


generated by the columns of B as K(B). Clearly, K(B) = {x € R"|x= » @;b;, ; 2 0}. 


i=] 
The cone we most frequently encounter is R° , the nonnegative orthant, where 
nxn 


R* = Kd) with I being the identity matrix in R 


The cone K induces an order relation on R" by x < y<y-xeK and 
x <y@y-xeK°® The dual, or adjoint, cone is K = {o € Ro, x) 20 Vx € K}, 
K 


which clearly satisfies the properties of a cone, and for a proper, sumplicial cone, 


(K') =K. Ifwe define Ko as K\{0}, then x e K° b,x) >O0VOE K,, and 


x € OK & (6, x) =0 for some @ € K,, where OK isthe boundary of K. 


A continuous function f; Dc R" > R’ is quasimonotone nondecreasing in x 
relative to the cone K if for x, y €¢ D, y— x € OK implies there existsa ¢ € Kj such 
that (6, y— x) =0 and , fly) — f{x)) = 0. This definition, from Elsner (1974), is fairly 
standard; however, for nonlinear functions some authors require x, y € K (see, for 
example, Heikkila, 1983). Fora linear function f{x) = Ax, A € R"”, the quasimonotone 
nondecreasing property reduces to: x € OK implies there existsa @ € Kj such that 
(d, x» =0 and <0, Ax) 2 0. 

An equivalent definition is f: R° > R" is quasimonotone nondecreasing if and only 
if f(x) is nondecreasing in x; forall 1j=1,...,n,1#J, so the following lemma 
characterizes quasimonotone nondecreasing linear operators. This result appears as an 
example in Lakshmikantham and Leela (1977b). 

Lemma A.1. A linear operator A € R®™ is quasimonotone nondecreasing 
relative to the nonnegative orthant if aj 20 for all 14). 

Proof. We prove this for n= 2. A similar argument works for n> 2. Let 
K = R2 = {(x), X2)|x1 2 0, x. 20}. Then K =K, since (, x) =0 implies 


1X1 + 2X2 2 0 for all x), x2 = 0, which in turn implies 0), 6220. If f(x) = Ax, where 


43, 4y2] . : : : 
A= is quasimonotone nondecreasing relative to K, then for x;=0 or 
do} 422 


X2 = 0, there existsa @ € Kj such that (b,x)=0 and (6, Ax) = 0. Now 
(), X) = 1X1 + 2X2 = 0 and (, Ax) = O)(a11X1 + a12X2) + O2(a21X1 + az2X2) 2 0. If x; = 0, 


X72 + 0, then 2X2 =0 so d2 =0. This implies, since 01 Se QO, that )1412X2 = 0, SO aj22 0. 


Similarly, az; 2 0, so the linear functions which are quasimonotone nondecreasing relative 
to K are precisely those matrices with nonnegative off-diagonal entries. 0 

We use the following well-known results. 

Lemma A.2. Given a nonsingular matrix B € R" andthe cone K(B), then 


K((B™)') = KG). 


Proof. Let x € K((B')') and y € K(B) so that x= = aid; and y= > Bibi, 
i=] 


i=] 


with a;, B; 20. Then (x, y)= y » O70 bia »S iB; > 0, so K((B")') c K(B). 
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Similarly, let x ¢ K((B™')'). Since the 6; are independent, x = y yi0;, but there 


i=1 

exists at least one yj < 0. For y=b; € K(B), it follows that (x, y)=y,;<0,so x¢ K(B). 
Hence, K(B) < K((B’')'), proving the lemma. 0 

Lemma A.3. Given a nonsingular matrix B € R" andthe cone K(B), then 
K(B) CR® if and only if K(B")>R°. 

Proof. Let K(B)c R* andlet ace R‘{ bearbitrary. Then Ba=x eR! since 
xX 1S a positive linear combination of the generators of K(B). Because B is invertible, it 
follows that a= Bx. As a € R® is arbitrary and expressible as a nonnegative linear 
combination of the columns of B™’, we conclude that K(B™') > R°. 


Conversely, assume K(B')>R®° and let a ¢ R®. Then there exists an 


x € R® such that B"'x=a,so x=Ba. Then B is nonnegative, otherwise it has some 
column b; which is not nonnegative, and letting a=[0),.. ., 0;-1, 1), Oj1,.. ., On]’ shows 
x ¢ R*, which isa contradiction. This proves the lemma. 0 

Since K(B) c R®° implies B is a nonnegative matrix, and since B’ is also 
nonnegative, then K((B™)') D> R* aswell. This gives a nice characterization of the 
inverses of nonnegative matrices. 

When we say that a matrix A € R"™” (ora linear operator A: R"— R’) is 
quasimonotone nondecreasing with respect to a cone K(B), we mean that C = B™'AB is 
quasimonotone nondecreasing (with respect to R’ ) even though A may not be. In 
Chapter III we show this result as a theorem of Heikkila. 

Since quasimonotonicity is related to nonnegative matrices, we next present some 


background on nonnegative matrices. 


B. NONNEGATIVE MATRICES 

The study of nonnegative matrices began with Perron’s (1907) presentation of a 
theory of positive matrices. Frobenius (1908, 1909, 1912) immediately extended this 
theory to nonnegative matrices, and the general theory has carried their names ever since. 
Here we present only their basic results required for our application. 

A positive matrix A € R®™" has a;>0 fori.j=1,...,n. A nonnegative matrix 
has a; 20. Amatnx A is essentially nonnegative (positive) if and only if aj =0 (> 0) 


for all 14+j. The following lemma justifies the term “essentially” nonnegative. 


Lemma B.1. Ifa matrix A € R"” is essentially nonnegative, then there exists a 
number r>0 such that A =A+trI is nonnegative. Furthermore, the spectrum 
o(A ) = o(A) +1, and the eigenvectors of A and A are equal. 


Proof. The first statement is obvious. Let Ax = Ax. Then Ax= (A + rI)x = 
Ax + rx = Ax + rx = (A + 1)x, completing the proof. 0 

Hence, we can shift the diagonal (and the spectrum) of an essentially nonnegative 
(positive) matrix to make it nonnegative (positive). 


A matrix A is reducible if it is permutation sunilar to a matrix 


A A 
= i. rm where Aj; and Az are square (and of positive dimension). 
22 


Otherwise, A is irreducible. A positive matrix is a special case of a nonnegative, 
irreducible matrix. The following theorem summarizes the basic results we need from the 
theory of nonnegative matrices. 

Theorem B.2. (Perron, 1907; Frobenius, 1912). A positive matrix A> 0 has an 
eigenvalue equal to the spectral radius of A. Associated with this eigenvalue is a positive 
eigenvector. Furthermore, A has no other nonnegative eigenvectors. The same 1s true 
for a nonnegative, reducible matrix A = 0. 

A nonnegative matrix A= 0 has a (real) eigenvalue equal to the spectral radius of 
A. Associated with this eigenvalue is a nonnegative eigenvector. 0 

For a proof of the above results see, for example, Gantmacher (1959), Varga 
(1962), or Hom and Johnson (1991). Since the transpose of a matrix with any of the 


above properties retains that particular property, the same conclusion can be made about 


10 


the left eigenvectors of A. Note that when comparmg two matrices asin A > B, or two 
vectors asin x > y, we refer to component-wise majorization. 

From Lemma B.1 and Theorem B.2 it follows that an essentially positive matrix 
has a real eigenvalue with greatest real part, with positive left and nght eigenvectors 
associated with this eigenvalue, and no other nonnegative eigenvectors. Similar 
conclusions follow for essentially nonnegative and irreducible essentially nonnegative 
matrices. 

Ifa lmear operator A: R"— R" is quasimonotone nondecreasing, then the matrix 
A € R" is essentially nonnegative. Therefore, if we ask when a square, linear operator 


A on R" is quasimonotone nondecreasing with respect to a cone in R‘., we are asking 


nxn 


the equivalent question of when the matrix A € R"™ is essentially nonnegative under a 
nonnegative change of basis. The next chapter introduces the motivating application, the 


theory of cone-valued vector Lyapunov functions. 
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IH. CONE-VALUED VECTOR LYAPUNOV FUNCTIONS 


The motivating problem for determining the quasimonotonicity of a square matrix 
with respect to a nonnegative cone in R" comes from stability theory in differential 
equations. Specifically, we seek such a cone when using vector Lyapunov functions to 
determine the stability of an equilibrium in a dynamical system. This chapter presents this 
problem for lmear comparison systems. 

We begin with a brief introduction to the theory of differential equations and 
stability. In order to develop the technique of vector Lyapunov functions, we also require 
results from the theory of differential inequalities. We present only the results from these 
fields needed to prove the major theorems of vector Lyapunov functions. In the case of 
well-known results, or those which will not be used directly, we present the results 


without proof. 


A. DIFFERENTIAL EQUATIONS , INEQUALITIES, AND STABILITY 
We consider a system of first-order differential equations 
x’ = f(t, x), x(to) = Xo, (A.1) 
where x’ = dx/dt, x € R", f Ix D—R™ where ICR and DCR” are open sets 
contaming to and Xo respectively, and f has certain continuity requirements. A 
solution x(t) to equation A.1 is a differentiable function x: J C R > R” satisfying the 
differential equation dx/dt = f{t, x) and the initial condition x(to) = Xo. We also consider 


the autonomous system 


x’ = f(x), x(to) = Xo, (A.2) 
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where the vector field is independent of time. 

Picard (1890) and Lindel6f (1894) showed that if f{t, x) is continuous on a closed 
set containing (to, Xo), and is uniformly Lipschitz continuous with respect to x on this 
set, then equation A.1 has a unique solution x(t) on some set [to, tp + a] for a> 0. 
Peano (1890) showed that without the assumption of Lipschitz continuity, then equation 
A.1 has at least one solution x(t) on t € [to, to + a]. Furthermore, a solution on an open 
subset of R™”’ containing (to, Xo) can be extended to the boundary of that set. Proofs of 
these results can be found, for example, in Hartman (1982), and these existence and 
uniqueness results are the basis for the theory of differential equations. 

To prove stability results about equilibria in equation A.1, we also need results 
from the theory of differential mequalities. We extend the concept of the standard 
derivative used in equation A.1] to the Dini derivative, defined as 


D_x(t) = Lum int (1/h)[x(t + h) — x(t)]. This is one of four Dini derivatives, and we define 
—0- 


the others as needed. We begin with scalar mequalities (m= 1 im equation A.1). The 
following results are found, for example, in Lakshmikantham and Leela (1969a) or Walter 
(1970). 

Theorem A.1. Let E be an open (t, x)-set in R* and fe C[E, R]. Assume 
that v, w € C[[to, tt + a), R] for some a > 0, and (t, v(t)), (t, w(t)) € E for 
t € [to, to + a). If v(to)< w(to) and if D_v(t) < f(t, v(t)) and D_wi(t) > f{t, w(t)) for 
t € (to, tp + a), then for t € [to, tp + a), v(t) < w(t). O 

For the scalar differential equation A.] (m= 1), we define the maximal solution on 


[to, to + a) asthe unique r(t) such that for every solution x(t) defined on [to, to + a), 
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x(t) <r(t) for t € [to, to + a). The minimal solution is defined analogously, and the 
following results presented for maximal solutions are also valid for minimal solutions, with 
appropriate changes. Under the hypothesis of Peano’s theorem, a maximal solution exists 
on the interval [to, to + a], and as before, if one exists on an open set, it can also be 
extended to the boundary of that set. 

It is useful when employing vector Lyapunov functions to compare a solution to a 
maximal solution over an interval. The followmg comparison theorem is basic to this idea. 

Theorem A.2. Let f © C[E, R], where E is an open (t, x)-set in R’, and let 
[to, to + a) be the largest interval in which the maximal solution r(t) of equation A.1 
exists. Let p € C[(to, to+ a), R], (t, p(t)) € E for t € [to, to + a), p(to) S Xo, and fora 
fixed Dini derivative, Dp(t) < f(t, p(t)) on t € [to, t) + a) except on possibly a countable 
subset of this interval. Then for t € [to, tp + a), p(t) <r(t). OC 

The question of the existence of a maximal solution for a system of differential 
inequalities was solved by Wazewski (1950). He extended the above result to the case 
m> 1 in equation A.1, and showed that a sufficient condition for a maximal solution to 
exist to equation A.1 is that f(x) be nondecreasing in x; for each t € [to, tp + a], a 
property which we defined in Chapter II as quasimonotone nondecreasing. This property 
was first recognized as being important to differential inequalities by Miller (1926) (see 
Walter, 1970). There are analogous definitions for creasing, decreasing, etc., as well as 
for mixed quasimonotone properties. 

Burton and Whyburn (1952) used the mixed quasimonotone property to prove the 


existence of what are known as minimax solutions to differential equations. These are 
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solutions which, for example, are minimal in the first k components and maximal in the 
last m—k components. The existence of minimax solutions leads to a family of 
comparison theorems for differential inequalities which can be used to prove results about 
stability of equilibria using vector Lyapunov functions. These include results about 
stability, instability, and conditional stability. We present only the comparison theorem for 
stability, sce we use stability for the motivating example of cone-valued vector 
Lyapunov functions. 

Theorem A.3. Let E be an open (t, x)-setin R™ , with fe C[E, R™). 
Suppose that f is quasimonotone nondecreasing in x, and let [to, to + a) be the largest 
interval of existence of the maximal solution r(t) of equation A.1. Let 
p € C[[to, to + a), R"], (t, p(t)) € E for t € [to, to + aw), and on this interval, for the Dini 
derivative, D_p(t) < f{t, p(t)). Then p(to) < xo mplhes p(t) <1(t) on this interval. 0 

This theorem applies for any Dini Derivative. The proof of this result, along with 
all the above results on differential inequalities, are in Lakshmikantham and Leela (1969a). 

We now present some basic definitions and concepts from stability theory which 
allow us to develop the technique of vector Lyapunov functions in its most simple setting. 

An equilibrium solution to equation A.1 is apomt € such that f(t, €)=0 for all 
t. Equation A.1 admits the trivial solution &=0 if f(t, 0)=0 for allt, and we frequently 
use this as our equilibrium as we can generally send an equilibrium point to the trivial 
solution via a change of coordinates. 

We wish to determine the stability of the equilibria of equation A.1, and we begin 


with some basic definitions of stability. The trivial solution of equation A. 1 is stable if for 
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every neighborhood U, of the origin and every tp 2 0, there is a neighborhood U2 ofthe 
origin such that xo € U2 implies x(t) € U; for all t 2 to, where x(t) is a solution 
satisfying x(to) = Xo. The trivial solution is asymptotically stable if it is stable and there 
exists a neighborhood U; ofthe origin such that xo € U3 implies that 

x(t) > 0 as to. 

We may express these definitions m terms of norms as follows. The trivial solution 
is stable if the solution x(t) with x(to) = xo is such that for every ¢ > 0, to 2 0, there 
exists a 6>0 such that ||xo||<6 implies ||x(t)||<e for all t >to. Similarly, the trivial 
solution is asymptotically stable if there exists a 6 >0 such that ||xo|| <6 implies that 
x(t) > 0 as to. 

The trivial solution in a linear system (when f(t, x) in equation A.1 is a linear 
function of x, ie., f(t, x) = A(x) where A € R™) is asymptotically stable when the 
eigenvalues of the matrix A all have negative real part for t 2 to. In this case we refer to 
the origin as a sink. In a nonlinear system, if the eigenvalues of the linearized system in a 
neighborhood of the equilibrium all have negattve real part, then that equilibrium is a sink. 

If the origin is a sink in a linear system, then for all solutions x(t), ||x(t)|| isa 
strictly decreasing function of t, where ||-|| 1s the Euclidean norm. Similarly) it is a well- 
known result (see, e.g., Hirsch and Smale, 1974; or Amold, 1989) that when the origin 1s 
asymptotically stable in a nonlinear system, there exists some norm such that ||x(t)|| 1s a 
decreasing function of t for all solutions x(t) startmg sufficiently near the ongin. 

These basic definitions of stability have numerous refinements. The two major 


classes of stability definitions are Lyapunov and Poisson, and among the logical statements 


ly 


allowed in the definitions are 17,017,969 possible types of stability (see Bushaw, 1969). 
Many of these conditions are meaningless, and Massera (1949, 1956) discusses those 
which are most often studied in terms of Lyapunov functions. Habets and Peiffer (1973) 
extend these to a classification which encompasses all possible types of stability. 

The above definitions of stability can be generalized for nonautonomous systems 
(for example, when 6 depends on to) and they can be extended to conditional stability 
(instability) definitions for stable (unstable) manifolds through the point (to, 0). There is a 
well-developed theory of vector Lyapunov functions for conditional stability which uses 
the minimax solutions of Burton and Whyburn, but we discuss only the most basic case. 
Furthermore, these stability results are local, and we do not present the conditions under 
which they are global. 

The definitions of stability and asymptotic stability used above rely on the norm of 
a solution being either bounded im any neighborhood of the trivial solution, or decreasing 
uniformly as t increases. Particularly for nonlinear systems, these conditions may be 
difficult to verify. Lyapunov recognized this, and in his 1892 dissertation he suggested 
that a function other than a norm could be used to determine stability (Lyapunov, 1907). 
Of course we now refer to such functions as Lyapunov functions, and while they are 
frequently difficult to find, they offer another powerful tool for classifying equilibria. 

Before introducing Lyapunov functions, we define the derivative of a function 
along a trajectory, in this case where the trajectory is a solution to the initial value 
problem. Fora solution x(t) to equation A.1 and a real-valued function V(t, x), we 


define (using the terminology of Wiggins, 1996) the orbital derivative of V along x as 
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V (t, x) = limsup (1/h)[V(t + b, x + hf{t, x)) — V(t, x)]. We now present a version of 
h—0* 


Lyapunov’s theorem from Hirsch and Smale (1974) for autonomous systems. A function 
V satisfying the hypothesis of the theorem is called a Lyapunov function. 
Theorem A.4. (Lyapunov, 1907). Let the autonomous system (equation A.2) 
x’(t) = f(x) admit the trivial solution f(0)=0. Let V(x) be a continuous function defined 
on a neighborhood U ofthe ongin such that 0 €e UCD. Let V be differentiable on 
U \{0}, with 
V(0O) = 0, and V(x)>0 if x #0. (A.3) 


Then 


V <0 on U\{0} (A.4) 


unplies the ongin is stable. Furthermore, 


V <0 on U\{0} (A.5) 

woplies the origin is asymptotically stable. 

Proof. Let U be any neighborhood of the ongin, and choose 56 such that 
B;(0) c U, where B;(0)is a 6-ball about the ongin. Let a be the minimum value of V 
on the boundary of B;(0), with a >0 by equation A.3. Let U; = {x € B;(0)|V(x) < a}. 
Then no solution starting in U, can meet the boundary of B;(0) since V is 
nonincreasing by equation A.4. Therefore, the ongin is stable. 

If equation A.5 holds as well then V is strictly decreasing on orbits in U\{0}. 
Let x(t) startin U,\{0} and suppose x(t,) > Z € Bs(0) for t, > 0. (By the 


compactness of Bs(0) such a sequence exists.) We now show that z=0. By the 


Ie 


continuity of V and equation A.5, V(x(t,)) — V(z) and V(x(t)) > V(z) for all t>0. If 
Z #0, let z(t) be a solution starting at z. For any s>0, V(z(s)) < V(zo), so for y(s) 
starting sufficiently near Zo, V(y(s)) < V(z). If we let y(0) = x(t,) for sufficiently large 
n, then V(x(t, + s)) < V(zZ), which is a contradiction. Therefore, z = 0 is the only 
possible limit point of {x(t)|t > 0}, and such a limit exists by the compactness of B;(0). 
This proves the theorem. OU 

While Lyapunov functions provide a useful technique for classifying equilibria, 
there is no known algorithm for findmg them. In mechanical systems or electrical circuits, 
potential or total energy is frequently a candidate, but in general no technique works 
uniformly. Even after 100 years, finding a Lyapunov function is still more of an art than a 
science. The following example is from Hirsch and Smale (1974). 


Example A.5. Consider the system 


x(t) = 2yz — 2y 
y(t) = 2+ x 
z'(t) =—z’. 


The origin is clearly an equilibrium, but if we linearize the system near the ongin via the 


Om —2 0 
Jacobian matnx then Jlp=|1  O 0}, which has eigenvalues 0 and +/2i. Weare 
0 O O 


unable to classify the origin via linearization, so we seek a Lyapunov function of the form 


V(x, y, Z)= ax’ + by’ + cz. Then 
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V = 2(axx’ + byy’ + czz’) = 2[2axy(z — 1) — bxy(z—1) — cz’]. Since we want V =O, 


letting c=1 and 2a=b gives V =-z’<0,s0 V=x’+2y'+Zz isa Lyapunov 


function, and the ongin is stable. However, the Lyapunov function is not strict since 


V = 0 on the xy-plane, and we cannot determine (with this choice of V) if the origin is 
asymptotically stable. 0 

While no known algorithm exists to find Lyapunov functions, they have proven 
very useful over the past one hundred years in determining the stability of equilibria. See, 
for example, Amold (1989), Bailey (1966), Bhatia and Szeg6 (1970), La Salle and 
Lefschetz (1961), Lehnigk (1966), and Yoshizawa (1966). 

In 1962, Bellman recognized that the requirement that a Lyapunov function be 
real-valued was too restrictive, and he proposed that a vector-valued function might 
provide more flexibility. He presented his results for square, linear systems and special 
cases of nonlmear systems. The following lemma is a special case of Theorem A.3 for 
linear systems, but we present it in its entirety as it provides necessary and sufficient 
conditions for the solution of a differential equation to be majorized by a solution of a 
differential inequality. In this lemma, A = [aij] € R"". The lemma is from Beckenbach 
and Bellman (1965). 

Lemma A.6. Let the system of differential equations 

dx/dt = Ax, x(0) = Xo 
have a solution x(t). Let y(t) satisfy the differential mequality dy/dt 2 Ay, y(0) = Xo. 


Then y(t)2x(t) for t20 ifand only if aj; 20 for 14). 
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Proof. The system dy/dt = Ax + g(t), y(0) =x» has a solution of the form 


t 
y(t) = xoe“! + | e'-9(s) ds, so we need to know when the elements of e* are 
0 


nonnegative fort>0. Since e“'=J+ At+A’t’?/2! + A°t’/3!+..., for small positive t 
e“‘xI+ At so that y(t) > x(t) implies a; > 0. 

Conversely, assume aj > 0 (the case aj 20 follows via a limiting procedure) and 
let e“t=(e*™)". For a fixed t, e”Y =1+ A(t/N) + A°(t/N)’/2! +... is positive (in the 
sense that all elements are positive for N sufficiently large). Since the product of positive 
matrices is positive, e“ is positive if a; > 0 for i#j. This proves the lemma. 0 

The condition that a; 20 for i#j is precisely the condition of f(x) = Ax being 
quasimonotone nondecreasing in Theorem A.3 for a lear function. Also, since a linear 
function is uniformly Lipschitz continuous, Picard’s theorem applies and the maximal and 
minimal solution are the unique solution x(t) = xoe“. 

We now State a stability theorem for vector Lyapunov functions, where we 
determine the stability of the trivial equilibrium solution to equation A. 1. 

Theorem A.7. (Lakshmikantham, 1965). Let g e C[R- x R!., R'], g(t, 0) =0 
for all t, let g(t, u) be quasimonotone nondecreasing in u for each t € R+, and let 


u’(t) = g(t, u), u(to) = Uo 2 0. (A.6) 
Let V € C{[R.x D, R{] and V(t, x) be locally Lipschitz in x, with Vi(t, x) > 0 as 
1=] 


|x|] + 0 for each t; and let f € C[R.x D, R"], ft, 0) = 0, and for (t, x) € R+ x D, 


V(t, x) < a(t, V(t, x). 


pup 


Then the stability of the trivial solution u=0 of equation A.6 implies the stability of the 
trivial solution x = 0 of equation A.1l. 0 
In the more general theorem of Lakshmikantham (1965), various types of stability 
and conditional stability are proven as well. Instead of presenting the proof here, we 
present instead an example from Lakshmikantham (1974) that demonstrates the usefulness 
of a vector Lyapunov function over that of a scalar Lyapunov function. 
Example A.8. Consider the system 
dx/dt = e“x + ysint — (x° + xy’)sin’t, (A.7a) 
dy/dt = xsint + e“y — (x’y + y’)sin’t. (A.7b) 
We would like to determine the stability of the trivial solution to this system, so we 


attempt to find a scalar Lyapunov function, trying V(x, y) = x° + y’. The best bound we 


can achieve with this choice of V is V < 2(e + |sint|)V. Clearly the trivial solution of 
du/dt = 2(e~ + |sint/)u is not stable, so we cannot conclude anything about the stability of 


the trivial solution of equation A.7 with this particular choice of scalar Lyapunov function. 


V(x, «i avi ary 


V ee 7 ane as a vector 
2 ? 


However, if we choose V(x, y) = 


Lyapunov function, then V 1 < 2(e~ + sint)V}), and V 2 <2(e' — sint)V2. Since the trivial 
solution u = 0, w= 0 to the system 

du/dt = 2(e~ + sint)u 

dw/dt = 2(e — sint)w 


is clearly stable, the trivial solution x = 0, y=0 to the system A.7 is also stable. O 


While we used the quantity »y Vi(t, x) as a measure, we could have used some 


i=] 


other measure such as max Vi(t, x) or Q(Vi(t, x), ..., Va(t, x)) where Q: R$ +R. 


and Q(u) is monotone nondecreasing in Uu. 

Vector Lyapunov functions offer a lot of flexibility when they can be found, but 
their true versatility is apparent in the case of conditional stability; 1.e., when the stable 
manifold of the equilibrium has dimension less than m. Here we present the most basic 
results of this theory m terms of conditional stability, which are from Lakshmikantham 
(1965). For completeness we review some of the terminology. 

We wish to determine the stability of the (generally trivial) equilibrium solution to 
equation A. 1, 

x’(t) = f{t, x), x(to) = Xo, to 2 0, (A.1) 
where x € R" and f{t, x) is defined and continuous in Rs x R". The trivial solution 
x = 0, where f(t, 0)=0 for all t 2 0, is conditionally equistable (where 5 im our 
definition of stability depends on ¢ and t,) if and only if there exists a manifold Ma-x 
through the ongin of dimension (m — k) such that for each ¢ > 0 and each tp 2 0 there 
exists a positive function 6 = O(to, €) that is contmuous in to for each ¢, where ||x(to)/| < 
& and x(to) € Mam+y imply ||x(t)||<« forall t 2 to. 

We consider the comparison system 


r’(t) = g(t, r), r(to) =m 20 (A.8) 
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where g(t, r) is quastmonotone nondecreasing in r € R". When g has this property, the 
solution r(t; to, ro) 1S maximal in the sense of component-wise majorization. We consider 
solutions of the form 
r(t; to, 1) =O Gi=1,2,...,k) (A.9a) 
T(t; to, fo) > O G=ktl,...,n) (A.9b) 
where m-—k corresponds to the dimension of the manifold Mn-t. 
The comparison system A.8 is equistable (where 6 in our definition of stability 


depends on € and to) if and only if for all ¢ > 0, to = 0, there exists a positive function 


6 = d(to, €), continuous m to for each £, such that sy r, <6 implies 


i=k+] 


n 


2 r(t; to, fo) <€ for t >to. Other stability and conditional stability properties are 


i=l 
formulated similarly. 

The vector Lyapunov function V(t, x) € C[R+ x R", R‘{] is locally Lipschitz 
continuous in x, and its orbital derivative along the trajectories x(t) of equation A. 1 


satisfies 
V(t, x) < a(t, V(t, x)) (A.10) 
where x(t) is any solution such that V(to, xo) < ro. By Theorem A.3 this ensures that 
V(t, x(t)) S r(t; to, 10) for t = to. 
The stable manifold of the origm, Mgt, 1s defined by the set of points for which 


Vit, x)=0 for 1=1,2,...,k<m (A.11) 
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The bound for the measure, b(r) is a continuous function, nondecreasing in r, 


such that 


b(r)>0 for r>0, and b(||x\|) < 3 Vi(t, x). (A.12) 


i=] 
We might also have the bound satisfy b(t) © as r — o, which is required for certain 
types of stability. 
While we have used the component-wise sum of the Lyapunov function as the 
measure, there are several equivalent choices we could have selected. In terms of this 


measure, we may require, for example, 


> Vit, x) > 0 as ||x|| > 0 foreach t20. (A. 13) 


i=l 
We now state the theorem on vector Lyapunov functions and conditional stability. 
Theorem A.9. (Lakshmikantham, 1965). Let assumptions A.10 through A.13 
hold. If the solution A.9 to equation A.8 is conditionally equistable, then so is the 
equilibrium solution of equation A. 1 


Proof. Let ¢>0. If ||x|| =, then from assumption A.12, we have 
b(e)< >) Vi(t, x). (A.14) 
i=l 


If the stability property holds, given b(e)>0 and to 2 0, there exists a positive function 


} = d(to, €), continuous in to for each ¢€, such that 


» r,(t; to, To) < b(€) for t= to, (A. 15) 


i=l 


26 


provided » r, <0. Let x(t) be any solution of equation A.1. Then it follows from 


1=k +1 


assumption A.10 that 


Vi(t, x(t)) < >. T(t; to, fo) for t = to (A.16) 
i=] i=] 
whenever 
De (te) a eae (A.17) 
i=] i=] 


Now choose r, (i= 1, 2,...,n) to satisfy 


ag ae for i= 1,2,...,k and (A.18) 
DS Eee. (A.19) 
i=k+1 


By equation A.17 and smce V(t, x)20 for 1=1,2,...,n, equation A.18 implies that 
x(to) € Mam-+ because of assumption A.11. From the monotonic property of b(r), 


assumption A.12, and equations A.17 and A.19, we deduce that ||x(to)|| < b’(8) = 8:. 


By assumption A.13, there exists a 52 = 62(to, 5) suchthat sup 3. Vi(to, x(to)) < 6. 
Ix(to ISS j=} 


Let 53 = min{6,, 52}. It then follows from the choices of i and 63 that x(to) € Mat 
and ||x(to)|| < 63 implies every solution x(t) satisfies equation A.16. Suppose, if possible, 
that a solution x(t) of equation A.1 satisfying that x(to) € May and ||x(to)|| < 53 is such 


that ||x(t)||=« for some t >to. Then from assumption A.12 and equations A.14, A.15, 


and A.16 follows the contradiction b(e)< pa V(t, x(t)) S be I(T; to, To) < b(e), 


i=] i=] 


2), 


which proves the theorem. 0 

Vector Lyapunov functions can also be used to determine certain boundedness 
properties of solutions of equation A.1, and boundedness theorems analogous to the 
above theorems for stability are given in Lakshmikantham (1965). We now give an 
example from Lakshmikantham (1965) of the usefulness of a vector Lyapunov function in 
determining the conditional stability of an equilibrim. 

Example A.10. Consider the system of differential equations 

dx/dt = (1 + cos(t))x + (1 — cos(t))y + (cos(t) — 1)z 

dy/dt = (1 —e“)x +(1+ e")y+(e" - 1)z 

dz/dt = (cos(t) — e“)x + (e" — cos(t))y + (e~ + cos(t))z. 
We select as the Lyapunov function 


V, (t,x, y,Z) Gey): 
Vit, x, y, Z)=| V,(t,x,y,z)| =| (x-ytz)’ 


3 
Since ~~ Viexty +27 +(x-y)y +(y—z) +(z—xy, we may use as the positive, 
1=1 


monotone function b(x, y, z) the square of the Euclidean norm, since 


3 
b(x, y,z)=x ty +2 < ss Vi(t, x, y, z), so condition A.12 is satisfied. Since 


i=1 
V i(t, x, y, Z) = 2(x + y — z)(dx/dt + dy/dt — dz/dt) = 4(x ty - zy 
V (t, X, y, Z) = 2(x — y + z)(dx/dt — dy/dt + dz/dt) = 4cos(t)(x — y + Z) 


V x(t, x, y, 2) = 2(-x + y + z)(—dx/dt + dy/dt + dz/dt) = 4e(—x + y +2)’, 
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the comparison system is defined by 
g(t, Vi, V2, V3) = 4Vi(t, x, y, Z) 
g(t, Vi, V2, V3) = 4cos(t)V2(t, x, y, Z) 
g3(t, Vi, V2, V3) = 4e V3(t, x, y, Z). 


The comparison system is therefore 


r, (t) ©. (Lesreeste) Ar, 
r, (t)| =r(t) = g(t, 1) =| g,(t,1,,6,1,)] = | 4cos(t)r, 
ee (t) £;(t,1,,0 .1,) oe 


Since each gj; is only a function of 1;, the function g is trivially quasimonotone 


nondecreasing, and satisfies V< g(t, V). 

Since the trivial solution to r,’(t) = g;(t, r) is clearly unstable, we must choose 
k = 1, so the imitial condition is r, = 0, and 1(t; to, 0) = 0. This forces the condition 
V(t, x, y, Z) = 0, which defines the two-dimensional manifold Mg-;) by the condition 
(x+y—z) =0 to be the pinne x +y=z. The solution to the comparison system is 

1(t) = 0, 

r(t) =r, exp[4(sin(t) — sin(to))], 


eS eet AG ac 3); 


3 
If ¢ > 0, then for 5 = min{e/(2e’), e/(2exp(4e °° )}, r,, +1, <6 implies >) 1(t)<e 


i=l 
for t >to. Since the conditional equistability condition is satisfied, the trivial solution to 


the original system is equistable if the initial point is on the manifold x +y=z. 0 


pe) 


The technique of vector Lyapunov functions is a proven technique in stability 
analysis in many fields, particularly for large scale dynamical systems and interconnected 
systems. See, for example, Michel and Miller (1977), Siljak (1978), Lakshmikantham, 


Matrosov, and Sivasundaram (1991), and Vidyasagar (1993). 


B. CONE-VALUED VECTOR LYAPUNOV FUNCTIONS 

While vector Lyapunov functions provide a great amount of flexibility over scalar 
Lyapunov functions, Lakshmikantham (1974) noted that the quasimonotonicity of the 
comparison system is not a necessary condition for the system to be stable. In particular, 
for a linear system, a matrix can still be a stability matrix (have all eigenvalues with real 
part less than zero) without having all the off-diagonal elements nonnegative. 

By requiring the vector Lyapunov function to have each V; 2 0, we have restricted 
it to the cone R{. Lakshmikantham and Leela (1977a) investigated the possibility of 
selecting a cone other than R‘ to overcome this limitation. The following theorems 
from Lakshmikantham and Leela (1977b) extend some of the previous results to cones 
other than R‘. The first is the comparison principle corresponding to Theorem A.3 
through the cone K. 

Theorem B.1. (Lakshmikantham and Leela, 1977b). Let f ¢ C[R. x R", R*] be 
quasimonotone nondecreasing in x relative to K foreach t € R.-, and let [to, ©), to 2 0, 
be the largest interval of existence for the maximal solution r(t; to, Xo) of equation A. 1 


relative to K. Further, let p € C[R+, R"] and D_p(t) < f(t, p(t)) for t2to. Then 
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p(to) S Xo implies p(t) : I(t; to, Xo) for t2to. O 


The quasimonotonicity of f(t, x) in x relative to the cone P does not necessarily 
imply the quasimonotonicity of f(t, x) in x relative to the cone Q when 
P <Q. However, if P <Q, then the order relations relative to P do imply the same 
order relations relattve to Q. We prove these observations in Section D, and from them 
comes the following corollary. 

Corollary B.2. Let P and Q betwoconesin R" such that Pc Q. Let the 


assumptions of Theorem B.2 hold with K replaced by P. Then p(to) < Xo implies 
p(t) < I(t; to, Xo) for t2to. O 


We now state the comparison theorems and stability results for cone-valued vector 
Lyapunov functions from Lakshmikantham and Leela (1977b). We begin with the system 
of differential equations 

x’(t) = f(t, x), x(to) = Xo, (A.1) 
where f € C[R. x D, R"]. If K isa cone in R", n< m, and the cone-valued vector 
Lyapunov function V € C[R. x D, K], we define for (t, x) € Rs x D, the orbital Dini 


derivative as D’ V(t, x) = limsup (1/h)[V(t + h, x + hf{t, x)) — V(t, x)]. The first result 
h-—>0* 


follows. 
Theorem B.3. (Lakshmikantham and Leela, 1977b). Assume that V(t, x) 
satisfies a local Lipschitz condition in x relative to K and for (t, x) € R: x D, 


D'V(t, x) s g(t, V(t, x)), (B.1) 


where g € C[R: x K, R"], and g(t, u) is quasimonotone in u with respect to K for 
each te R.. If r(t; to, uo) is the maximal solution of u’(t) = g(t, u), u(O) = uo relative to 


K, and x(t; to, Xo) 1s any solution of equation A.1 such that V(to, xo) & Uo, then on the 


common interval of existence 
V(t, x(t; to, Xo)) ss T(t; to, Uo). (B.2) 
Proof. Let x(t) = x(t; to, xo) be any Petition as above. Set p(t) = V(t, x(t)). For 
small h > 0, since V(t, x) is locally Lipschitz in x with respect to K, then 
p(t + h) — p(t) S Li}x(t + h) — x(t) — bit, x(t))]] + V(t + h, x(t) + hf(t, x(t))) — V(t, x(t). 
From this and equation B.1 follows the inequality D"p(t) = g(t, p(t)). Applying Theorem 


B.1 gives the conclusion B.2. 0 

The following variant of this theorem offers more flexibility in applications, and its 
proof follows from Corollary B.2. 

Theorem B.4. (Lakshmikantham and Leela, 1977b). Let P and Q be two cones 
in R" such that P< Q. Suppose that V € C[R: x D, Q], V(t, x) satisfies a local 


Lipschitz condition in x relativeto P, and D’V(t, x) < g(t, Vit, x)) for (, x) € R. x D, 


where g € C[R: x D, R"] and g(t, u) is quasimonotone nondecreasing in u relative to 
P for each t € R.. If r(t; to, uo) is the maximal solution as in Theorem B.3 relatrveto P 


and x(t; to, Xo) 1s any solution of equation A.1 such that V(to, xo) < Uo, then 


V(t, x(t; to, Xo)) S r(t; to, Xo) (B.3) 


EY: 


on the common interval of existence of r(t; to, Uo) and x(t; to, Uo). When Q=R°%, 
relation B.3 implies the componentwise inequalities V(t, x(t; to, Xo)) <r(t; to, Uo). CO 

We now State the first stability results for cone-valued vector Lyapunov functions. 
The first theorem extends the result of Theorem A.7 for equistability. Its proof follows 
from that of Theorem A.9 and Theorem B.3. 

Theorem B.5. (Lakshmikantham and Leela, 1977b). Let the assumptions of 
Theorem B.3 hold. Let f{t,0)=0 and g(t, 0) =0, and assume that for some 
do € Ky and (t, x) € R: x D, 

b(IIx||) < (do, Vit, x)) S a(t, |]xl), (B.4) 

where a € C[R. x [0, p), R+], b € C[[0, p), R-+], a(t, 0) = 0, b(0) =0, a(t, u) and b(u) are 
increasing in u, and p = sup{||x||/x € D}. Let the trivial solution u=0 of u’(t) = g(t, u) 
be o-equistable, that is, given ¢ > 0, to 2 0, there exists a 5 = O(to, €) > O such that 
(do, Uo» < 5 implies (do, r(t; to, Uo)) <€& for t >to. Then the trivial solution x=0 of 
equation A.1 is equistable. 0 

The following theorem again increases our flexibility when employing cone-valued 
vector Lyapunov functions. The proof follows from those of Theorems A.9 and B.4. 

Theorem B.6. (Lakshmikantham and Leela, 1977b). Let the assumptions of 
Theorem B.4 hold, and let f{t,0)=0 and g(t,0)=0. Assume equation B.4 1s satisfied 


for some oo € Q, and that the trivial solution u =0 of u’(t) = g(t, u) is o-equistable, 


with 9 € Qj. Then the trivial solution x =0 of equation A.1 is equistable. 0 


U2 
2 


it KR then oo ere, 1]', and Theorem B.5 reduces to Theorem A.7. 


In this case, condition B.4 becomes b((|x||)< >) Vi(t, x) S a(t, ||x/|). We could again use 


i=l 
measures other than the component-wise sum of the vector Lyapunov function. In 
Theorem B.6, if Pc Q= RR‘, we can remove the requirement for quasimonotonicity if 
we Select an appropriate cone P depending on the nature of g(t, u) as shown in the 
followmg example from Lakshmikantham and Leela (1977b). 
Example B.7. Consider the system 
u)'(t) = aya; + apU2 = g)(t, Uy, U2), Ui(to) = Ud), (B.5a) 
U2’(t) = a2)U) + axU2 = 22(t, U1, U2), Ueto) = U2. (B.5b) 
If Q=R*, and az, and aj2 are not nonnegative, then the function g(t, u) is not 
quasimonotone nondecreasing in u = (uj, U2) relative to Q. Hence, the differential 
inequalities 
D'Vi(t, x) S$ gi(t, Vit, x), V2(t, x)), (B.6a) 
D'Vxt, x) S$ g(t, Vi(t, x), V2(t, x)) (B.6b) 
do not yield the componentwise estimates of V(t, x(t)) im terms of the solution of 
equation B.5. However, if there exist two numbers a and B suchthat 0<B<«a and 
oraz, + Oar > ay; + a1, (B.7a) 
B’a21 + Baz = Ban + ap, (B.7b) 
then these conditions hold with no restriction on the nonnegativity of a2; and aj2. Let 


the cone PC Q=R* be defined by P= {u € R* |Buy Su; < ap}. The boundaries for 


this cone are Gu2=u, and Bu,=u). On the first boundary, let 6 =(-—1/a, 1) so 
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((-I/a, 1), (ui, ui/a)) = 0, and ((-I/a, 1), (ayy; + ay2uj/a@, a2, + a22u;/a~)) > 0 for all 
u #0. This reduces to condition B.7a, and we similarly obtain condition B.7b. Thus, if 
the inequalities B.6 are relative to P, we obtain component-wise estimates on V as 

Vi(t, x(t)) S r(t; to, V(to, xo)) (B.8) 
by Theorem B.4. If ai2, a2; 2 0, then equation B.8 is the one we would obtain through the 
standard method of vector Lyapunov functions. 0 

While it appears to be a nontrivial exercise to construct an appropriate cone for a 
system of differential equations, this method still carries much merit as it provides a further 
increase 1n our flexibility to determine the stability of equilibria. 

These are the most basic results from the theory of cone-valued vector Lyapunov 
functions we require to motivate the problem of finding a nonnegative cone with respect 
to which a linear differential operator on R" is quasimonotone nondecreasing. Further 
directions in the theory of vector Lyapunov functions have been explored recently. In 
particular, using higher derivatives of vector Lyapunov functions, K6ksal and 
Lakshmikantham (1996) showed how to find a particular cone with respect to which a 
given comparison system (resulting from taking higher derivatives of a decrescent 
Lyapunov function) is quasimonotone. 

Furthermore, there has been much research on nonlinear comparison systems. See, 
for example, Hatvani (1984), Deimling and Lakshmikantham (1990), Lakshmikantham, 
Leela, and Ram Mohan Rao (1991), and Lakshmikantham and Papageorgiou (1994). 
However, we focus on finding cones for vector Lyapunov functions when the comparison 


system is linear. 


C. THE QUASIMONOTONICITY OF LINEAR DIFFERENTIAL SYSTEMS 

The problem of finding a cone in R" with respect to which a given linear operator 
is quasimonotone nondecreasing (or essentially nonnegative) and many related problems 
have been addressed in various forms in the literature in recent years. For example, 
Vandergraft (1968) following Birkoff (1967) gave sufficient conditions for such a cone to 
exist without further restrictions on the cone K. However, Heikkila first addressed this 
problem for the application of vector Lyapunov functions, with the requirement that the 
cone be proper, simplicial, and nonnegative. He stated this to be sufficient for the 
application of vector Lyapunov functions, and m the next section, we justify this claim. In 
this section, we summarize Heikkila’s results. 


nxn 


Consider a lmear mapping A € R°” andacone K generated by n vectors b; m 
R" such that 

K = {wyb, +... + Wab,|w; 2 O}. (C.1) 
For C.1 to define a proper, simplicial cone in R", it is necessary and sufficient that the b; 
be linearly independent. The first result gives necessary and sufficient conditions for a 
linear operator to be quasimonotone with respect to a cone in R". 

Theorem C.1. (Heikkild, 1980). A linear mapping A: R" > R" is 
quasimonotone nondecreasing relative to the cone generated by a basis {bi,..., ba} of 
R" if and only if the matrix of A relative to this basis has all off-diagonal elements 
nonnegative. 


Proof. Let A =[a,;] and C =[cj] be the nxn matrices of A relative to the 


standard basis of R" and the basis {b,,..., b,} respectively. Then C =B ‘AB, where 
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B < R™ has b; asits jth column. Let $? be the ith row of B’', so 
(i, bj) = 04, (C.2) 
and 
(Qi, Ab;) =~ Cij. (63) 


Equations C.2 and C.3 imply that 6; € K; for each i=1,...,n. 


Assume ¢;; 20 for i#j, and let u € OK be given. Thenif u= >) wyb,, at least 


jel 


one coefficient w; must be zero. Then (0, u) = » W;<;, bj) = 0 and 
jel 


(oi, Au) = 5° wj(gi, Abj)= >> wjej 20 from equations C.2 and C.3. Thus, the 


. . 


jel j=l 
condition cj 20 for 1+j implies the quasimonotonicity of A relative to K. 

Conversely, let at least one off-diagonal element of C, say cj, is negative. Then 
by equation C.3, (;, Ab;) <0. Since the mapping u — (@;, Au) is continuous, there 
exists a 6>0 such that ||u —b;||<6 implies 


{o;, Au) < 0. (C.4) 


n 


In particular, equation C.4 holds for u € OK given by u= ys wyb;, with w; = 0, 


k=1 
w; = 1, and for k #1, j, Wk = 0/(al|b, ||). Since all the coefficients wx, except wi, are 
positive, then for any b € K; for which <b, u) = 0, it follows that <b, b,) =0 for all 


k #1. Thus, b must be of the form b=a9; for a> 0, so by equation C.4, 


(b, Au) = ag), Au) < 0. This shows A is not quasimonotone nondecreasing relative to 
K unless cj 20 for 14#j, completing the proof. 0 

Salzmann (1972) proved the same result for nonnegative matrices, and he 
presented it m terms of positive operators on simplicial cones. Furthermore, in his 
discussion of matrices with invariant cones, Vandergraft (1968) proved a similar result. An 
immediate consequence is the following. 

Corollary C.2. (Heikkila, 1980). Ifa lmear mappmg A: R° > R* has n 
Imearly ndependent eigenvectors b;,...b,, then A is quasimonotone nondecreasing 
relative to the cone generated by these eigenvectors. 

Proof. In this case, the matnx of A relattveto {b;,...,b,} is diagonal. OD 

We can also express a necessary condition for monotonicity in terms of the 
eigenvectors of A. This result follows from the result from Perron (1907) and Frobenius 
(1912) on the theory of nonnegative matrices, which states that a nonnegative matrix has a 
nonnegative eigenvector corresponding to a real eigenvalue of greatest modulus (see 
Theorem II.B.2). 

Corollary C.3. (Heikkila, 1980). Ifa lear mappmg A: R° > R" is 
quasimonotone nondecreasing relative to a cone of type C.1, then A has a non-zero 
eigenvector in this cone, and hence a real eigenvalue. 

Proof. If A, = BAB is quasimonotone nondecreasmg where B = [b;,.. ., Dal, 
choose a positive number r so large that C=A,+rI isnonnegative. Then C hasa 
nonnegative eigenvector x; 20,so A, (and hence A) has an eigenvector 


yi Bx, € K(B). a 
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We now consider linear differential systems. By the quasimonotonicity of the 
linear system 
u(t) = A(t)u, t 2 0, (C5) 
relative to a cone K we mean that foreach t 20 the lmear mapping A(t) is 
quasimonotone nondecreasing relative to K. The next result for cones contained in R ° 
follows from Corollanies C.2 and C.3. 
Corollary C.4. (Heikkilé, 1980). For the quasimonotonicity of the system C.5 


relative to some cone in R‘ it is necessary that for each t 20 at least one eigenvector 
of A(t) belongsto R‘{, and sufficient that A(t) has n lmearly independent eigenvectors 


in R‘°, which do not depend on t. O 


nxn 


Given a nonsingular matrix B € R°””’, the mapping u = Bv transforms equation 
C.5 into the equivalent form 
v(t) =BA(t)Bv, t 2 0. (C.6) 

But B”'A(t)B is the matrix of the linear mapping A(t) in the basis of R" formed by the 
column vectors of B. Moreover, if B is nonnegative, these vectors belong to R°. 
Hence, the next result follows from Theorem C.1. 

Theorem C.5. (Heikkilé, 1980). There exists a cone K in R‘ generated by n 
linearly independent vectors, relative to which the system C.5 is quasimonotone 
dondeereastie if and only if there exists a nonsingular matrix B such that the off-diagonal 


elements of the coefficient matrix of the system C.6 are nonnegative. Moreover, the 


column vectors of B generate K. U 


This special case of considermg conesin R‘ lets us apply Theorem B.6 when the 


vector Lyapunov function satisfies a differential equality. 


D. STABILITY THROUGH LINEAR COMPARISON SYSTEMS 

The results of Heikkilé (1980) presented in Section C address the problem of 
finding a cone with respect to which a linear comparison system is quasimonotone, and 
provide a necessary condition and a sufficient condition for.such a cone to exist. Kéksal 
and Lakshmikantham (1996) constructed cones for comparison systems generated by 
taking derivatives of Lyapunov functions, and Kéksal and Fausett (1995) extended these 
results using generalized eigenvectors. 

Heikkilé (1980) states (without proof) that for stability results it is sufficient to 
consider Bae cones. This section gives our proof of this result, and shows it is 
sufficient to consider square, linear comparison systems. 

We wish to determine the stability properties of the origin as a solution to 

x’(t) = f(t, x) (A.1) 
where f € C[Rs x R™ > R"] and f{t, 0)=0. To do this we seek a vector Lyapunov 
function V(t, x) € C[R+ x R", K] which is locally Lipschitz in x with respect to a cone 


K c R’, and which satisfies the differential inequality 


V < a(t, V) (D.1) 


where V is the orbital derivative of V along trajectories x(t). If g(t, u) is 


quasimonotone nondecreasing with respect to K, the stability of the trivial solution of 
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u(t) = g(t, u) (D.2) 
may be used to determine the stability of the origin x(t)=0 in equation A.1. We require 


certain conditions on the Lyapunov function V for this result to hold, for example, 


n 


b(Ix|l)< >) [Vil S a({lxll) (D.3) 


i=] 
where a and b are asin Theorem B.5. 

We let the comparison system A.3 be a square, linear system, or g(t, u) = A(t)u 
where A: R: x R" > R’, or for all t, Ae R™. Let Ve R® so the quasimonotonicity 
of A(t) implies its off-diagonal elements are nonnegative, and the mequality D.1 is 
component-wise majorization (otherwise we can change bases so V € R° J. 

If A(t) is a stability matrix, but does not possess the required quasimonotonicity 
properties, we seek a cone K CR" other than R‘ for which inequality D.1 holds, 
relative to which A(t) 1s quasimonotone nondecreasing, and which contains V(t, x) 
while maintaining its locally Lipschitz property. The following lemma discusses the effect 
of a lmear transformation on an inequality. 


Lemma D.1. If x < y, then for a nonsingular matrix B € RR", Bx & | ‘By, 
1 


where K(B) c R” is the cone generated by the columns of B. 


Proof. Since Bx < By <= By ~— Bx € K(B) and (y—x) € R°, then from the 


K(B 
linearity of B, B(y — x) € K(B) = {u= >) ajbjla; = 0, b; is a column of B}. Letting 


i=] 


ai =(y — x); = 0 completes the proof. 0 


4] 


We can transform a differential system v’(t) = Av, with a nonsingular matrix B, 
where y = Bv, so y’(t)=BAB"'y. Such a transformation is called a Lyapunov 
Transformation and the stability of the origin in one system implies the stability of the 
origin in the other. In this case, we say the systems are equivalent in the sense of 
Lyapunov (see Gantmacher, 1959). 

We now show that a simple coordinate transformation cannot be used to construct 
a cone which satisfies the hypothesis of Theorem B.6 when the original system does not 
possess the required quasimonotonicity property. 


Proposition D.2. Assume a vector Lyapunov function satisfies the inequality 
V < AW, (D.4) 
but that A(t) is not quasimonotone nondecreasing with respect to R{. Then there does 
not exist a Lyapunov transformation Y= B™'V for which Y aah C(t)Y and C(t) is 


quasimonotone nondecreasing with respect to K(B’’). 

Proof. If such a transformation existed, then with C(t) = B™'A(t)B, multiplying 
inequality D.4 by B™ and applying Lemma D.1 gives 

Y=B'V <s_ B'A(t)V=B BC(t)B V = C(t)Y. 
K(B) 

However, in order that C(t) be quasimonotone nondecreasing with respect to K(B™'), by 
Theorem C. 1 it is necessary and sufficient that BC(t)B™' have its off-diagonal elements 
nonnegative. But BC(t)B™' = A(t), which, by assumption, is not quasimonotone 


nondecreasing with respect to R?. C 
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Hence, the problem is to find an appropriate cone with respect to which the 
function A(t) is quasimonotone nondecreasing, and for which the differential inequality 
D.1 and all other required properties are preserved. The following case, where the 


inequality D.1 is really an equality, is frequently found in applications. 


If the orbital derivative of a Vector Lyapunov function V € R° satisfies V =AV 
where A is a stability matrix, then the following results give conditions under which the 
hypotheses of Theorem B.6 are satisfied, and the solution x(t)=0 to equation A.1 isa 


stable equilibrium. 
First, it is evident that if V = AV then B V s BAV for all nonsingular 


transformations B and all cones K. We also need the following result. 


Lemma D.3. If V(t, x) € R{ is locally Lipschitz in x with respect to R‘, then 


nxn 


for any nonsingular matrix B € R°”, Y = BV is locally Lipschitz in x with respect to the 
cone K(B). 

Proof. Evidently, V(t, x) € R implies Y = BV e€ K(B). Since V 1s locally 
Lipschitz with respect to R‘, this implies there exists an L € R‘ such that 
Lx - y|| - (V(t, x) - V(t, y)) € R®. Letting L =BL e€ K(B) yields | 
L\Ix - yll - (Y(t, x) - Y(t y)) € KB). 0 

We also require the result that the Lipschitz property is preserved through cone 
containment. 


Lemma D.4. If V(t, x) € P is locally Lipschitz in x with respect to the cone 


PCR‘, and PCQ where QC R" isa cone, then V(t, x) is locally Lipschitz in x with 
respect to Q. 

Proof. If V(t, x) is locally Lipschitz in x with respect to P implies there exists 
an Le P such that L||x — y| -— (V(t, x) — V(t, y)) € P, then L € Q and 
L||x — yl - (V(t, x) - Vit, y)) € Q. 0 

The following result appears as a remark in Lakshmikantham and Leela (1977). 


Lemma D.5. If P and Q areconesin R" with Pc Q, then x S y implies 


Proof. x Sy > y-xeP > y-xEQ =x Sy. 0 


We now demonstrate the well-known result that if P and Q are conesin R" with 
P <Q, that the quasimonotonicity of a mapping A € R™" with respect to Q does not 


imply the quasimonotonicity of A with respect to P. Let P=R*% and let Q bethe 


| 1 0 0 1 
cone generated by the columns of B = . : , hence P< Q. Let A= ; z 


-—3 1 
which is not quasimonotone nondecreasing with respect to P. Since BAB = ; } 


then by Theorem C.] A is quastmonotone nondecreasing with respect to Q. 
Similarly, if P c Q are cones in R’, the quasimonotonicity of a mapping with 


respect to P does not imply the quasimonotonicity of the mapping with respect to Q. 


—2 
Lettnmg A= ; | and choosing P and Q as above, A is quasimonotone 


S. 
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=—§ 2 
nondecreasing with respect to P, but B-'AB = SB A , so by Theorem C.1, A is not 


quasimonotone nondecreasing with respect to Q. 
Following Heikkild (1980) and Kéksal and Fausett (1995) we present the 
following result. 


Proposition D.6. Suppose a Vector Lyapunov function V(t, x) € R" satisfies 


V = A(t)V where A(t) is a stability matrix which is not quasimonotone nondecreasing 
with respect to R{. Ifthe generalized left eigenvectors of A(t) are all contained in R° , 
there exists a transformation Y=B™V and acone K=R” such that the hypotheses of 


Theorem B.6 are satisfied. 

Proof. The matrix A has as its Jordan canonical form J = B”'AB, where the rows 
of B™ are the generalized left eigenvectors of A and the columns of B are the 
generalized eigenvectors of A (see Hom and Johnson, 1991). By assumption, B™ isa 


nonnegative matrix, so K(B™') c R®. Hence, the transformation Y =B™'V yields the 


differential equation Y =JY. Then Y € K(B”) is locally Lipschitz in x with respect to 


K(B™') by Lemma D.4, and with respect to R" by Lemma D.5. Since J is 


quasimonotone nondecreasing with respect to R', and Y < JY, weuse Y asthe 
R 


vector Lyapunov function with K = R“, satisfying the conditions of Theorem B.7. The 


only remaining detail is the measure used to bound Y, so the following remark completes 


the proof. 0 
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Remark D.7. Since P< Q implies Q’ c P’, we cannot, in general, find a 


bE Q, =R\{0} such that (>, BV) preserves the original measure 
(bo, Vy= >° V; with oo =[1,..., 1]’. This is because $ = Boo is not generally in 
i=l 


Q* =R"\{0}, since the nonnegativity of B™’ does not imply the nonnegativity of B. 


However, had we used the equivalent measure max|Vi| = max Vi, then using as the new 
1 1 
measure max |Yjj, since each Yj; is a positive lear combination of the V;, preserves the 
1 


result of Theorem B.6 by scalmg a and b appropriately in equation B.4. 0 

While we used the Jordan form to achieve quasimonotonicity with respect to R°, 
this result may be extended to any case where S = B ‘AB, S has nonnegative off-diagonal 
elements, and the columns of B™' are contained in R“. This follows from Theorem C.5 
of Heikkilé (1980), who further shows that this construction is always possible in the case 
n=2 if A has one eigenvector in the first quadrant and at most one of a;2 and az; 1s 
negative (in this case it is necessary and sufficient that both eigenvectors be nonnegative). 
This is an application of Corollary C.5 to square systems with differential equalities. 


In the case where we do not begin with an equality in equation, D.1, this 


construction is generally not possible. Instead, given V < A(t)V, we now must find a 
R 


cone with respect to which A is quasimonotone nondecreasing, on which V is defined 


and locally Lipschitz in x, and on which the inequality D.1 still holds. The result of 


K6ksal and Lakshmikantham (1996) is a special case of this. 
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We define the following four properties for a linear comparison system for a vector 
Lyapunov function V andacone K. 
(P1) V is defined on K; 


(P2) V is locally Lipschitz with respect to K; 
(P3) V < A(t)V; and 
(P4) A(t) is quasimonotone nondecreasing with respect to K. 

If we find a cone K =P for which properties (P1) through (P3) hold, then in 
order to find a cone K =Q for which properties (P1) through (P4) hold it is necessary 
that P< Q. This is because ifa cone Q containing P with respect to which A(t) is 
quasimonotone nondecreasing exists, then by Lemmas D.4 and D.5, properties (P1) 


through (P3) hold for Q, and we can apply Theorem B.6. 


This is essentially what was done mm Proposition D.2 where the vector Lyapunov 


function Y € K(B) C R°. satisfied Y - | JY, and J was quasimonotone 


nondecreasing with respect to R!. Lettmg K(B)=P and R‘ =Q satisfies properties 
(P1) through (P4). 

We first state an obvious result for finding such a cone, which may be useful m 
applications. The proof follows from the above remarks that P c Q. 

Proposition D.8. If V € P satisfies properties (P1) through (P3) for K = P, and 
ifthe cone K(B) generated by the generalized eigenvectors of A contains P, then 


properties (P1) through (P4) hold for K(B). O 
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This generalizes to any transformation S = B™'AB where S. has nonnegative off- 
diagonal elements and K(B) contams P. Following Heikkilé (1980) we next show that 
to find a cone containmg R‘ with respect to which A is quasimonotone nondecreasing, 
it is necessary and sufficient to find a cone contained in R* with respect to which A’ is 
quasimonotone nondecreasing. We use Lemma II.A.3 which states that K(B) CR‘ if 
and only if K(B"') > R®. Since it is presumably easier to find a cone contained in R° 
(the matrix of its generators is nonnegative) than one containing R‘{ (as we require) this 
proposition is useful in applications. 

Proposition D.9. The matnx A is quasimonotone nondecreasing with respect to 
acone K(B")>R° ifand only if A’ is quasimonotone nondecreasing with respect to a 
cone K(B')cR°. 

Proof. Let A’ be quasimonotone nondecreasing with respect to a cone 
K(B') CR®. Then by Theorem C.1, C =(B')'A'B' has all of its off-diagonal elements 
nonnegative, as does C' = BAB”. Since K(B')c R® implies K(B) CR®, then 
K(B')>R® and A is quasimonotone nondecreasing with respect to K(B’’). 


Reversing the argument completes the proof. C 
Hence we can apply Corollary C.4 and Theorem C.5 of Heikkila (1980) to the 
search for an appropriate cone, as the generalized eigenvectors of A’, or the generalized 


left eigenvectors of A, bemg contained m R® isa sufficient condition for the existence 


of a cone for which property (P4) holds. 
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Because it is the case most frequently encountered in applications, we have 
assumed P= R_ to establish properties (P1) through (P4). Evidently, if PCR? the 
above results still hold; however, if P>R‘ then the construction is not as simple. 

Since for two cones P and Q,P CQ implies Q' cP, thenif P ¢ R®, instead 
of seeking a set of vectors bj € R° such that (B')'A'B' has nonnegative off-diagonal 
elements, where the b; are the columns of B’, we must find vectors b; € P’ to construct 
an appropriate cone. We therefore conclude the discussion of square linear comparison 
systems with the following extension of Proposition D.9. 

Theorem D.10. Given a vector Lyapunov function V and acone PCR” for 
which properties (P1) through (P3) hold, in order to find a cone QDP with respect to 
which the matrix A is quasimonotone nondecreasing it is necessary and sufficient to find 
m independent vectors b; € P such that if B' =[b, ... ba], then C =(B')'A'B' has 
nonnegative off-diagonal elements. 

Proof. If b; ¢ P are such that C=(B')'A’B’ has nonnegative off-diagonal 
elements, then so does C' = BAB”, and A is quasimonotone nondecreasing with respect 
to K(B™') by Theorem C.1. Since K(B™') = K(B'),, and since K(B') CP’ implies 
K(B') >P, then K(B™')>P. Reversing the argument completes the proof. 0 


Proposition D.9, where P= R£ is a special case of the above result, since 


R° =R°.. Finding an invertible matrix whose columns are contained in P’ is perhaps 


+ 
less difficult than finding one whose columns generate a cone contaming P. The next 


result follows from Corollaries C.3 and C.4 of Heikkila (1980). 
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Corollary D.11. For property (P4) to hold relative to some cone Q > P, where P, 
Q, V, and A(t) are as above, it is necessary that at least one eigenvector of A’(t) belong 
to P and sufficient that A‘(t) has n eigenvectors in P’, which do not depend on t. 0 

To show quasimonotonicity with respect to an appropriate cone, we only show 
such a cone exists. In an equivalent approach we could use the cone to transform the 
vector Lyapunov function and the comparison system as we did in Proposition D.2. This 
technique now works, since the cone used to make the transformation is contained in our 
original cone, and the inequality is preserved. The following corollary is an extension of 
Theorem 2.8.3 of Lakshmikantham, Matrosov, and Sivasundaram (1991). 

Corollary D.12. Ifa vector Lyapunov function V exists such that properties 
(P1) through (P3) hold for some cone P, and if a matrix B exists as m Theorem D.10, 
then the vector Lyapunov function Y =(B')'V gives the same conclusion about the 
stability of the trivial solution to Equation A.1. 0 

This approach has the advantage of directly yielding the stable manifold of the 
origin in the case of conditional stability, but in the usual case it involves an additional 
computational step. Since the ideas are equivalent, we continue the approach of showing 
such a cone exists. 

We continue our analysis of linear comparison systems considering the case where 
the system is rectangular instead of square. Let a nonnegative vector Lyapunov function 
V(t, x) be such that 


V: R.xR"SR?, (D.5) 


whose orbital derivative satisfies inequality D.4, V < A(t)V. Let A(t) € R'™ with 
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5 
>) 


A= (D.6) 


=) 
© 


where A € RA e RP*4 0 isa q Xp zero matnix, and 0 isa q x q zero matrix, 


with p+q =n. 


Such a comparison system is possible when the orbital derivative satisfies V ; < 0 


for p<j<n. If we have equality in any of these components, or if V j= 0 for some j, 
this is a conserved quantity for our original system A.1, and the theory of conservative 
systems to applies to our analysis. While it 1s necessary that a conserved quantity be non- 
constant on open sets, a component of a non-trivial vector Lyapunov function certainly 
satisfies this requirement, so this more powerful theory applies to the problem. Hence, we 


assume the orbital derivative satisfies inequalities in its last q components. 


This rectangular p x n system can be treated as a square system D.4 where A is 
given by equation D.6. Since we address such systems in Chapter IV, where we give a 
complete theory for reducible, square matrices, there is no need to consider the smaller, 


rectangular system. 


We similarly dismiss the case where the matrix A in equation D.5 is of the form 


>> 
a 


5] 


where A € R®? and the last q columns of A contain zeroes. A is not a stability 
matrix since o(A) = o( A ) U {0;,..., Og}, and A may not possess quasimonotonicity 


properties, even if A does. However, lettng Y;= V; for 1<i<p yields the square 


system Y <A Y, to which we may apply the theory for square comparison systems, 

regardless of the nature of A. Hence, we gain nothing by keeping the last q 

components of V, so we may discard them, leaving the imbedded p x p square system. 
We have now established the setting for the problem of finding a cone with respect 

to which a linear comparison system is quasimonotone nondecreasing. It is sufficient to 

consider square linear systems, and considering the transpose of the comparison system, it 

is sufficient to find a cone contamed m the nonnegative orthant. The next two chapters 


present our solution to this problem. 
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IV. THE QUASIMONOTONICITY OF A SQUARE, LINEAR OPERATOR 
WITH RESPECT TO A NONNEGATIVE CONE: THE REAL SPECTRUM 


For a matrix A € R™" we address the problem “when is A quasimonotone 
nondecreasing with respect to a cone in R? ?” (or equivalently, when is it essentially 
nonnegative under a nonnegative change of basis?). In Chapter III we motivated this 
question with the application of determining stability in dynamical systems via the 
technique of cone-valued vector Lyapunov functions, and in Chapter VI we present 
further applications of this problem. 

Using Perron-Frobenms theory, Heikkild (1980) showed a necessary condition (the 
matrix A has a nonnegative first eigenvector) and a sufficient condition (the matrix A 
has all nonnegative eigenvectors) for such a cone to exist (see Corollary I.C.2) where the 
first eigenvector is the nonnegative eigenvector x, associated with the (real) eigenvalue 
A, of greatest real part. 

The next two chapters present our solution to this problem. Using constructive 
techniques, we bring together the necessary and sufficient conditions Heikkild gave for 
n> 2, and we address the complex spectrum for the first tume in this setting. This chapter 
addresses the case where the spectrum of A isreal, or o(A) CR, and Chapter V 


addresses the general spectrum. 


A. MATRICES WITH A POSITIVE FIRST EIGENVECTOR 


We begin with the case where the first eigenvector is positive, or x; > 0. We use a 


sequence of changes of basis An+1=B—~ AmBn with K(B,) C R°, which ensures the 


change of basis is always nonnegative. The following lemma is the key to the 


construction. 


Lemma A.1. Let B= [x, b2,..., bh] € R"™" be nonsingular with x € (R‘)° and 


ay 
1 


T 


B” = Zz 
a3 
0, 
Then there exists a unit basis vector e, € R" such that the matrix 


iE 
©, 


p2 |e: 
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has a nonnegative mverse D € R‘"", where (+) indicates an appropmiate choice of sign 
for each 9). 

Proof. Since x € (R‘)° and (6;,x)=0 for i=2,...,n, then +0; ¢ R°. 

Hence, each $; has components of both signs. Since each hyperplane $; contains 
x € (R°¢ )°, then each of these hyperplanes intersects at least (n — 1) of the coordiate 
hyperplanes e; in R‘. (This is evident from the mixed signs in 9; for i=2,...,n.) 
Therefore, there are (n — 1) hyperplanes which each intersect at least (n — 1) of the 
(e)); =e; AR*. By the pigeonhole principle, at least one (e;); intersects each 6; . 
Let this unit basis vector be e,. Since we require <(¢;, d;)= 1 where d; is a column of D 
for i=2,...,n, to ensure d; € (e,); we select the appropriate sign for each 9;. 


Therefore, D™ has as its inverse D = [di, do, ..., da], which is nonnegative. 0 
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We note that the kth rowof D is e,, and that d) = x/x,. 


We now construct acone in R' with respect to which a matnx A is 


quasimonotone nondecreasing. We begin with the case where A is diagonalizable, and 


the first eigenvalue A, is sumple. In this case we can construct the matrix 


T 
ef, 


a °: (A. 1a) 


6: 
where the b; are left eigenvectors of A and K((B™')')> R° by Lemma A.1. Then 


B =[x), bo, . . ., b,] is nonnegative, and 


dr, ae Ain 
| 2 0 
A, =B AB= (A.2) 
0 Xr. 
Ifany a,, <0 for j=2,...,n, then we use another change of basis of the form 
l Sn 
] 0 
B, = (A. 1b) 
0 ] 
where the &; are all nonnegative, and 
Ay Aip 7 Any 
ei d, 0 i 
[Me BEINGS (A.3) 
0 Xr 


sD 


Simce ay = a1; + &(A1 — Aj), and A, > Aj, we can make aj; nonnegative, so A> is 
quasimonotone nondecreasing. 

In the case where 2, is not a simple eigenvalue, this construction still works, 
because for Aj = Ai, 41; =0 in equation A.2. To show this, since 4,;=e, Ab; and b; is 
in the eigenspace of A; =A, with [bj], =0 (since b; is orthogonal to all of the left 
eigenvectors ; where Aj # Aj, andto e, )’ Hence 4);=e, Ab; = 0. 

When A is not diagonalizable, Koksal and Fausett (1996) extended Heikkila’s 
results using generalized eigenvectors. Following this, we use the “almost diagonal” form 
of A, or the canonical form of A with arbitrarily small ¢ > 0 on the super-diagonal (see, 
for example, Horn and Johnson, 1991). Then, for negative 41;, aij = 44; + &(A1 — Aj) or 
aij = a4 + &(A1 — Aj) — €§)-1, so that aj; 20 for arbitrarily small ¢. Since the only other 
nonzero off-diagonals of A, are ¢ > 0, A> is again quasimonotone. 

The results of the previous two paragraphs combine in the case where A, 1s the 
eigenvalue which makes A ieteiee We summarize our results as follows. 

Theorem A.2. Let A € R"™ have a real spectrum and let a; <0 for some 1#J. 
In order for A to be quasimonotone nondecreasing with respect to a Tae cone, it 
is necessary that the first eigenvector x, = 0, and sufficient that x, >0. OU 

The disadvantage of the above construction is that we must find all of the 
generalized left eigenvectors of the matrix A in to construct the matrix A, (and the cone 
K(BB,)). However, since we only need to know whether such a cone exists for stability 


results, it may be sufficient to compute only the first eigenvector. 
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Since in the case where x, has some negative component we know no cone exists, 


and when x; 1s positive a cone always exists, the only remaining case for a real spectrum 
is when x; € OR°. 


B. REDUCIBLE MATRICES WITH A NONNEGATIVE FIRST 
EIGENVECTOR 


This section addresses the case where x; € OR‘, or when x, 20 has one or 
more components equal to zero. We assume x, € (R* )° ®0*%, where 0° indicates the 
zero vector in R*, and p+q=n. Because x; € OR’ we cannot necessarily apply 


Lemma A.1 to the left eigenvectors of A to obtain a nonnegative change of basis. 
Since under a (suitably ordered) nonnegative change of basis B, with x, € K(B), 


xP Bux e R?’ @©0". From this and Theorem I1.B.2, it follows that if A; = BAB= 


A. A.|. er ee ae 
e ay is quasumonotone, then if A,, is irreducible, it is no larger than p x p. 
ZZ 


Furthermore, if A,, is irreducible, then it can be put into the upper triangular form of 


J 


_, |, Where B> has the form of equation 
ORs. 


equation A.3 via a transformation B™ = 


A.1. Since this transformation does not effect the nonnegativity of A ,,, in this case the 


quasimonotone matrix A; = B ‘AB has a reduced form with a q x p zero block. From 
this follows our first result about reducible matrices, but we first clarify the particular form 
the reducible matrix must have. 


For a reducible matrix A with a first eigenvector x; € OR®, when A is 


permuted so that x, has all ofits zeros mm its last q components and A hasa q xp 


a7, 


lower-left zero block, we define A as reduced. We treat matrices which are not reduced 


when x, is ordered as above with the irreducible case. 


Ai, Ai 


Theorem B.1. A reduced matrix A = 
0 A, 


| is quasimonotone with 


respect to a nonnegative cone only if A;,; and A. are quasimonotone nondecreasing 
with respect to nonnegative cones. 


Proof. The previous remarks address the case where the quasimonotone form 


q + + 
An Ai Aj; 


A ASeeAS | 
BAB = A,= es res is reduced. Consider the matrix A, =| 0 Aj, Aj, 
22 (=o VA 


where Aj are quasimonotone and irreducible, and A; 20. Let the diagonal blocks 


have dimensions r x r,s x s, andt xt respectively, where r<p,t<q,r+s>p, and 


s+t>q. Since the largest eigenvalue , of Az is either an eigenvalue of AZ, or 


Aj,, then A, has associated with this eigenvalue an eigenvector which is nonnegative in 
the last s+t components. Hence, A, has a nonnegative eigenvector associated with 
Ul;. The case where this eigenvector is positive is the case of A q above, and when A22 
is reduced, we can repeat this argument. Since the result of this theorem is trivial in the 
case p = 1 (which is the only case required to prove Theorem C.2 in the next section) 
then the case where A>, is not reduced follows from that theorem O 

We now construct a cone when A is reduced, with A), € R°*? and Az € R™. 
Assume that o(Ai,) = {Ai >A22...2A)} and o(Az) = {Wi >b22...2 Ug} where 


A, >, and A>. has a positive first eigenvector. 
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B 0 
First let Bo = \ where B,, and B22 have the form of equation A.1, so 
22 


Xr, Cee, 
0 Xr 
Ae 8. AB, = ° 
0 oe 
0 Ly 
has upper-left and lower-right blocks with the form of equation A.3. While this 


transformation does not leave Ajz unchanged, we can now let A. be arbitrary. 
I B ss 
Next let B, = S a” where B)> = " ; tag with y; 20 and 


A> = B;'A,B,. This only changes the first row of A,,, where 


Apt > Arpe + Yi(A1 — bi) = aipe1, and Aips) > Aipy — YiBj + ¥/(Ar — Wy) = aij for 
j=2,...,q. Smce A; > yj, aj; can be made arbitrarily large. 

] 0 
some 


2 


D,, 0 
Now let B> = ate where D,)= with 06; = 0, and 


=D, l 

A; = B;'A,B,. While Bz is not nonnegative, we show that the product BoB,Bz2 still is, 
but first we examine the effect of this transformation. The off-diagonals in the first row of 
A2 are unchanged under B32, as are the last q rows. For 1=2,...,p,j=1,...,4q, 


Aip+j FD aAiptj) + 0)a1;, and since aj; 1s arbitrarily large, a;,+; can be made nonnegative for 


arbitrarily small 6;. Furthermore, for 1=2,...,p3;j=2,...,p31¥#J, aj 5 d;0; > 0. 


a7 


Pp 
Finally, for 1=2,...,p;aiy bb 6fA1 -Ai- ener ], and smce A, >A, andthe o& can 


k=2 
be made arbitrarily small, then a;,;20 as well. Hence, A; is quasimonotone. 


: ‘ - B,,D,, B,,B,, 
We now show that B = BoB,B2 is nonnegative. Since B= : 


0 BS 


we only consider B,;D,;. If we permute A,; so that when constructing B;;, e, =e; in 


1 0 
_ | X2 Dee 
Lemma A.1, then B;; = : , Where bj; 20 and x;>0,simce Aj; hasa 
x by OF 


positive first eigenvector. Only the first column of B;; changes under B,,D;;, and 


Pp 
Gh xX - >5 ;b, 20 since x;>0 andthe 8; are arbitrarily small. Therefore, in this 
i 


case, a nonnegative cone always exists. 

We now show we can relax the assumption that A;; and A22 have first 
eigenvectors which are strictly positive. If A;,; and Az are each quasimonotone 
nondecreasing with respect to nonnegative cones, then by the above construction, each 
can be put in a quasimonotone form with arbitrarily small nonzero elements below the 
main diagonal. In this case, the above construction still works for the matrix A. We 


summarize our results. 


12 


A,, A 
Theorem B.2. For a reduced matrix A = - ‘A | with real spectrum and 


22 
first eigenvector x; € (R?’ )° © 0° to be quasimonotone nondecreasing with respect to a 


nonnegative cone, it is necessary that A;; and A22 are quasimonotone nondecreasing 
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with respect to nonnegative cones, and sufficient that A, > u;, where A, and wu, are the 
first eigenvalues of A;,; and A>. respectively. 0 

The only reduced matrices which we have not included in this theorem are those 
where [; = A,. If both Aj, and A». are in upper-triangular quasimonotone form with 


4, = Ay, this construction may or may not be possible, depending on the sub-matrix Aj3. 


r 
The following example with n= 2 illustrates this. When A> pu, A= F "6 
rT 


] 
is quasimonotone with respect to K(B) where B= F fl for sufficiently large y, 


regardless of the sign of a:2. Also, if a;2 20, then A is quasimonotone with respect to 
K(), regardless of 1 and pp. However, if u2A and aj2<0,then A isnot 
quasimonotone with respect to any nonnegative cone. (The case u >A violates the 
necessary condition of the first eigenvector being nonnegative.) Furthermore, if a), = 0, A 
has two nonnegative generalized eigenvectors, but if aj2 <0, A has a generalized 
eigenvector x2 ¢ R*. 

Using the same construction as we did when A, > L, we formulate the following 
sufficient condition for a cone to exist, which we do not describe in terms of generalized 


eigenvectors, but instead in terms of the sub-matrix Aj2. 


Ay Ay 


with real spectrum and 
ae. 


Corollary B.3. For a reduced matrix A= 


first eigenvector x; € (R°’ )° © 0" to be quasimonotone nondecreasing with respect to a 


nonnegative cone, it is necessary that A,, and A22 are each quasimonotone with respect 


6] 


to nonnegative cones. Furthermore, if Ai; and A» are in upper-tnangular 
quastmonotone form with aj, =A; and ap+15+1 = 1 = Ai, where 1; is a simple 
eigenvalue of Azz, then it is sufficient that aj,+.20 for i=1,...,p. 0 

We can formulate similar sufficient conditions from this construction if 1 is not a 
simple eigenvalue of Az. Given the non-uniqueness of generalized eigenvectors, we have 
not produced a condition in this case that is based on generalized eigenvectors, although in 
the case where n= 2 the sufficient condition in Corollary B.3 is also necessary, and the 


requirement for both generalized eigenvectors to be in R~ is evident. 


C. IRREDUCIBLE MATRICES WITH A NONNEGATIVE FIRST 
EIGENVECTOR 


This section addresses the case where x; € OR‘, but A 1s either irreducible or 


7 —2 -7 7 —2 -7 
not reduced. Consider the matrices A=|4 1 —-9] and A=|6 -1 —-9 , which 
-1 J 1 -l1 1 3 


have A; = 5, Az = 2.894, A3 = 1.106, and x, =[1, 1, Oo}. A 1s quasimonotone 


| G2 ets 
nondecreasing with respect to K(B) where B=/1 3. 17), yet A isnot 
0 05 | 


quasimonotone nondecreasing with respect to any nonnegative cone. 
From the discussion preceeding Theorem B.1, if A has a quastmonotone form via 


a nonnegative change of basis, then some such matrix A, is reduced with a lower-left zero 


An Ap 


block. Furthermore, if A, -| neon 
22 


and Aj); has a positive first eigenvector (for 
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example, if it is reducible) then we can put A), in the form of equation A.3 using a 

transformation of the form of equation A.1. If 4; 1s suople, we can use the techniques of 

Section B for reduced matrices to ensure Ajz is nonnegative. This proves the following. 
Lemma C.1. Ifa matnx A with real spectrum and first eigenvector 


x; € (R?.)° © 0% is quasimonotone nondecreasing with respect to a nonnegative cone, 


then it is quasimonotone nondecreasing with respect to a nonnegative cone for which x; 
is an extremal. 0 

We cannot generally construct a nonnegative cone with respect to which A is 
quasimonotone as we did in Sections A and B (since the example in this section shows 


that such a cone does not always exist) but we can reduce the dimension of the problem by 


1) 0 
one by deflating the matrix with the first eigenvector as follows. Let Bo = h : 
xy 
where I is the identity matrix of dimension n—- 1, and x, =[1, &,..., &, 0, .. .,0]’, so 


May] fA ay here 4. = dfor i>p, a. = 
aA | , where a, = aj — aijC;, and for 1>p, a, = aj. 


A, = B:'AB -| 
1 0 0 : 0 ie 


Since Theorem B. 1! 1s trivial in the case p = 1, A; 1s quasimonotone with respect 
to a nonnegative cone only if A. 1s, so this is a necessary condition for A to be 
quasimonotone with respect to a nonnegative cone. Furthermore, since A, 2 A; for 1= 2, 
...,N (1.e., Aj is an eigenvalue of A>) then using a transformation B, of the form of 
equation A.lb, we can make A, = B,'A,B, quasimonotone as well. 

The cone K(Bo) is a subcone of R? © R24, where we replaced one of the first p 


basis vectors with x,. The transformation we showed replaced e,, but the choice of 
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which basis vector to replace can be varied by permuting the first p components of x). It 
is only necessary that one such transformation produce a matrix Az which is 
quasimonotone with respect to a nonnegative cone in order for the above construction to 
produce a cone for the matrix A. We summarize these results as follows. 

Theorem C.2. Let A € R® have a real spectrum with a first eigenvector 


x; € (R® )° © 0% and let A be irreducible (or not reduced). Then for A to be 


quasimonotone nondecreasing with respect to a nonnegative cone it is necessary that for 


1) O 
= ¢ ) 7 for some permutation x, € (R°)° ®0" of x, and for 
xy 


Anan: , 

ho Selby = FA ” | that the matnix A., be quasimonotone nondecreasing with 
a2 

respect to a nonnegative cone. In this case, it is sufficient that A, be a simple eigenvalue 

of A. 0 


When A, is not simple the problem is the same as the one discussed in the reduced 


di? 


" | and the quasimonotonicity depends on the sign 


r 
case, where, for example, A = 


of aj. 

The choice of which permutation of x; to select is not obvious, but there are only 
p such choices that lead to distinct cones. One solution is to compute the second 
eigenvector x2 of A, associated with A2, which is the eigenpair required for A222 to be 
quasimonotone. Since the last n— 1 components of B~'x: are the first eigenvector of 
Ax», then finding a B™' such that B™'x 20 is necessary, and that the last n— 1 


components of B™'x,> 0 is sufficient, to determine that such a cone exists. 
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V. THE QUASIMONOTONICITY OF A SQUARE, LINEAR OPERATOR WITH 
RESPECT TO A NONNEGATIVE CONE: THE GENERAL SPECTRUM 
A. IRREDUCIBLE, ESSENTIALLY NONNEGATIVE MATRICES. 

In this chapter we discuss the problem of finding a cone with respect to which the 
matrix A is quasimonotone nondecreasing in the case where A has a general (complex) 
spectrum, or when o(A) is not strictly real. Because A is real, its complex eigenvalues 
and eigenvectors occur in conjugate pairs; however, unlike with the real spectrum, the 
real-diagonal form of A is not quasimonotone (it has diagonal blocks which necessarily 
have a negative off-diagonal element) so we can no longer use the diagonal form. 

Additionally, a matrix A with a general spectrum is not always sumilar to a 
quasimonotone matrix, much less through a nonnegative change of basis. The problem of 
determining when a given set of complex numbers is the spectrum of a nonnegative matrix 
is known as the nonnegative inverse eigenvalue problem, and we discuss it in the next 
section. Since a matrix with a real spectrum is always similar to a quasimonotone matrix 
(its diagonal, or Jordan canonical form) this was not an issue in Chapter IV, and the only 
task was to produce a nonnegative change of basis. However, in this chapter, the 
nonnegative inverse eigenvalue problem plays a significant role in our solution. 

We have some amount of flexibility over the nonnegative inverse eigenvalue 
problem in that the matrices we seek are essentially nonnegative, and we can make them 
nonnegative by shifting the diagonal (and the spectrum) by a positive quantity r. Since 


BAB is essentially nonnegative if and only if B-'AB + rl is nonnegative for sufficiently 
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large r, and since B'(A+r1)B = B ‘AB +11, for this problem to have a solution the 
nonnegative verse eigenvalue problem must have a solution for the set o(A) + r. 
We now show that the matrices we need to consider have an irreducible, 


A, Ay 


is essentially 
Om” 


essentially nonnegative form. Otherwise, if B-AB = A,= 


nonnegative, then the shifted spectra o(A,,) +r and o(A22)+r must both solve the 
nonnegative inverse eigenvalue problem. In this case we transform Aj using the 
construction in Chapter IV, although if A,,; and A22 are not reducible, we cannot take 
advantage of the upper-triangular structure they had in the case of the real spectrum. We 
use this, and the methods of Chapter IV, to discard a portion of the real part of o(A), and 


we consider only that part which is required to form a solution to the inverse eigenvalue 


where 


I 
problem. For example, if o(A22) C R, by using transformations B = : 
Ze 


B22 is formed using the techniques of Chapter IV, we only need the theory for the 
complex spectrum to treat Aj). 

Since the final quasimonotone form is wrreducible, it has a positive first 
eigenvector, so A having a positive first eigenvector is a necessary condition for this 
problem to have a solution. (Or asin A, above, x; must have at least enough positive 
components to correspond to the dimension of Aj.) Furthermore, from Frobenius 
(19 12) the first eigenvalue A, must be simple or the quasimonotone form is reducible. 

Since the (shifted) spectrum of A must solve the nonnegative inverse eigenvalue 


problem for A to be quasimonotone nondecreasing with respect to a nonnegative cone, 
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and it is necessary that the first eigenvector of A be positive when the quasimonotone 
form, Ag, 1s irreducible. We now show that these conditions are sufficient. 

Without loss of generality, assume A, is not cyclic (otherwise we can further shift 
its spectrum so that it is non-cyclic). Let D be the real-diagonal canonical form of both 
A and A,, soD=B,' AByp = B, A,B. Not only are the spectra of A and A, the same, 
but the matrices are similar, so that the eigenvalues have the same geometric multiplicity 
as well. Hence, when we require the shifted spectrum of A to solve the nonnegative 
inverse eigenvalue problem, this includes the geometric multiplicity of the eigenvalues. 
Here, Bo =[X1,..., Xn] and B,; =[bi,...,b,] where the x; and b; are the real parts of 
the (generalized) eigenvectors of A and A, respectively, and x; and b, are the 
positive first eigenvectors associated with the simple eigenvalue 2, of greatest real part. 

While A, = B,B,' ABoB ,’, it is not generally true that BoB;’ > 0, as required. 
However, we further change bases using A, as n->oo without changing A, to make 
the change of basis nonnegative. Let B = BoB,’ AZ, so 


A,=B ‘AB = (A, ) BB >" ABoB | A,- Since A, 1s irreducible, non-cyclic, and has a 


simple first eigenvalue, A, — [bi,...,b:] as n>. Therefore, as n > «, 


re | 
B —> BoB - [bi, ..., Di] = Bo ; | = [x},..., X:] > 0, so A is quasimonotone 


nondecreasing with respect to a nonnegative cone. 
Hence, given a positive first eigenvector, we have reduced the problem to the 


nonnegative inverse eigenvalue problem,. We summarize our result as follows. 
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Theorem A.1. For a matrix A € R°™ with complex eigenvalues to be similar to 
an irreducible, essentially nonnegative matrix by a nonnegative change of basis, it is 
necessary and sufficient that A have a positive first eigenvector (associated with the 
simple real eigenvalue of greatest real part) and that the shifted spectrum of A solve the 
nonnegative inverse eigenvalue problem for an irreducible matrix. OD 

In the case where the final quasimonatone form is reducible, Az. as well as Aj, 
may have a non-real spectrum, in which case we can apply the above theorem to each of 
them individually, and then use the techniques of Chapter ITV to make Aj2 nonnegative. In 
the case where o(Az22) is strictly real we can use the techniques of Chapter IV to “throw 
out” this real part of the spectrum, and construct Aj, using the above techniques. The 
following example shows, however, that we cannot simply “throw out” all but 4, from 
the real part of the spectrum. 


The matrix A = isa 4x4 circulant matrix (each row is a shift of 


—_- Oo OO O 
oo .o = 
SS) a — 
oO = Oo © 


the previous one; see Davis, 1979), which is a type of cyclic matrix (its eigenvalues are the 
four roots of unity). Since o(A) = {+], +1}, this set is a solution to the nonnegative 
inverse eigenvalue problem. However, we cannot discard the eigenvalue A, =-—1, since 
the set o = {1,1,-i} is not the spectrum of any nonnegative (or quasimonotone) matrix. 


Furthermore, the left eigenvector associated with A, has two components of each sign, so 
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J 0 
that the transformation required to throw out 4, B= 62 ) cannot have a 
4 


nonnegative inverse B. 
When a matrix has a reducible quasimonotone form under a nonnegative change of 
basis we cannot discard all but the essential real part of the spectrum prior to constructing 


A,. Using transformations of the form B™ = for some collection of n—-p 


left eigenvectors, requires B to be nonnegative, which is not true in general. However, in 
the cases where this construction can be done a priori, then Theorem A.1 gives sufficient 
conditions for our problem to have a solution. 

To complete this topic, we present a brief discussion of the nonnegative inverse 


eigenvalue problem. 


B. THE NONNEGATIVE INVERSE EIGENVALUE PROBLEM 

There are many types of inverse eigenvalue problems. Given a set of numbers and 
a structure, they ask when that set of numbers is the spectrum of a matrix with the given 
structure. For a thorough survey of inverse eigenvalue problems, see Chu (1998). The 
structure we require is essential nonnegativity, or nonnegativity for the shifted spectrum. 
The nonnegative inverse eigenvalue problem is one of the classic unsolved problems from 


the theory of linear algebra (see, for example, Horn and Johnson, 1991) although some 
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sufficient and some necessary conditions have been given for a set of complex numbers to 
be the spectrum of a nonnegative matrix. 

The nonnegative mverse eigenvalue problem dates back to Perron and Frobenius, 
and has been approached in numerous forms and variants since then. Necessary conditions 
for real spectra were given by Sulemmanova (1949). His results were improved by 
Karpelevich (1951), Perfect (1952, 1953, 1955), Suleimanova (1965), and Fiedler (1972). 
Fiedler (1974) extended these results with the additional requirement that the matrix be 
symmetric, and Friedland and Karlin (1975) solved a eauan of the problem known as the 
nonnegative additive inverse eigenvalue problem. Most of the results for real spectra are 
refinements of Suleimanova’s necessary condition that the absolute value of the sum of the 
negative eigenvalues may not exceed A; > 0. Results of this type are always obtainable 
for essentially nonnegative matrices since we are allowed to shift the spectrum, and as 
mentioned earlier, the diagonal matrix of any real spectrum is essentially nonnegative. 

For complex (not strictly real) spectra, the problem is more difficult. Because of 
the close relation to nonnegative matrices (see von Mises, 1931) it has been addressed in 
various forms for nonnegative stochastic matrices, starting with Fréchet (1933), 
Romanovsky (1936), Dmitriev and Dynkin (1945), Taussky (1948), and Brauer (1952). 
Brauer presents a summary of the previous results, m addition to his own results and 
methods of proof. Friedland and Melkman (1979) resolve the question in the case of 
nonnegative Jacobi matrices, and numerous bounds for the spectra of nonnegative 
matrices exist (see, for example, Ostrowski and Schneider, 1960) as well as results on the 


first eigenvector of nonnegative matrices (see, for example, Minc, 1970; or Ashley, 1987), 
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but these latter results require a priori knowledge of the matrix. Additionally, Friedland 
(1977) extended his previous results for real spectra to complex spectra for the additive 
and multiplicative nonnegative inverse eigenvalue problem. 

Of the more recent results on the nonnegative inverse eigenvalue problem for 
complex spectra, as it applies to our problem, we mention those which provide sufficient 
conditions that may be used by the analyst, for example, to draw a stability conclusion in 
the method of vector Lyapunov functions. Kellogg (1971) addressed the problem of when 
a set of complex numbers is the spectrum of an essentially positive matrix (hence, for our 


application, certainly an irreducible essentially nonnegative matrix). He found a sufficient 
condition to be Re (A;) — Re(Aj+:) > aS; Im(A;+1) when the eigenvalues are ordered by 
their real parts, where 1, is the first (real) eigenvalue, and the A; are the eigenvalues 
with positive imagimary parts for i> 1. The other real eigenvalues do not affect this 
condition. Certainly, the condition is not necessary, as the set { V2 +6, +i} is the 
spectrum of a positive matrix for all € > 0. 


Friedland (1978) addressed the problem for eventually nonnegative matrices 


(where A“ is nonnegative for all k > M for some M >1) and showed that if Ss peal 


i=] 
for some p € {1],2,...} and if A, is the only positive 4; for i=1,...,n, then A is 
eventually nonnegative. This supports the conjecture that if k= 1, then {A;} isthe 


spectrum of a nonnegative matrix, although no proof has been found yet. Furthermore, if 
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A 20, A‘ 20, it is necessary that > M > 0. Similarly, Lowey and London (1978) gave 


i=] 


the further necessary condition that oa a esp ll a for alien ey 2 ee 


i=l i=l 

Boyle, ez. al., (1991) characterized via symbolic dynamics the sets of complex 
numbers which are the nonzero portion of the spectrum of a nonnegative matrix; however, 
their result is conjecture. Their results are important in that they are presented for 
primitive matrices (those matrices A 20 for which A‘ >0 for some k). For the 
nonnegative inverse eigenvalue problem, this is sufficient as the generalization to 
irreducible or nonnegative matrices is easily determined. More significant, however, is 
that for our stability application, since a primitive matrix is a special case of an irreducible 
nonnegative matrix, this is all we need to apply the techniques of the previous section. 

Further results, for example the sufficient conditions of T. J. Laffey or those of 
Koltracht, Newman, and Xiao (1993) (where Boyle’s conjecture is shown to be true for 
n < 5), are discussed in Berman and Plemmons (1994). There are probably enough 
necessary and sufficient conditions in the literature that most spectra encountered can be 
determined to be solutions of the nonnegative inverse eigenvalue problem (or not), but as 


long as the problem is unsolved, the pathological case cannot be discounted. 


TZ 


VL FURTHER APPLICATIONS 


This chapter presents additional applications that can benefit from the results of 
this dissertation. 

Boyd, et. al., (1994) present the problem of positive orthant stabilizability, which 
is a special case of the “hit and hold” problem from control theory (see Berman, Neumann, 
and Stern, 1991). Given a linear differential system x’ = Ax where A € R"", when does 
x(0) € R‘ imply x(t)e R{ forall t20, and x(t) 0as t > 0? It is a necessary and 
sufficient condition that there exist a diagonal matrix P>0O such that PA’ + AP <0, and 
that A be essentially nonnegative. 

Ohta, Maeda, and Kodama (1984) discuss a similar problem where the state 
variable is not required to be nonnegative, in which case the problem of finding an 
essentially nonnegative matrix is the nonnegative inverse eigenvalue problem for the 
shifted spectrum. They also present the positive realization problem, which is similar to 
the problem of positive orthant stabilizability, and for which essentially nonnegative 
matrices play an important role. 

In the problem of positive orthant stabilizability, we cannot arbitranly shift the 
spectrum of a linear dynamic system without changing the problem. However, if the 
matrix A isnot essentially nonnegative, but its spectrum satisfies the necessary 
conditions, then our techniques can find a cone contained in the positive orthant which can 
be stabilized under the original system. Our techniques can also be used to find a cone 


containing the positive orthant which can be stabilized under the original system. 


73 


A related problem from control theory in Stern (1980) is that of asymptotic 
holdability, where a trajectory must be held arbitrarily close to a linear subset of R" (not 
necessarily the nonnegative orthant) for any finite time interval. This leads to positively 
invariant cones, which Stern (1982) characterized for cones other than R*. These are 
cones K such that when x(t) is the solution to x’ = Ax, then x(0) € K implies x(t) e K 
for all t > 0. A sufficient condition in this case is the essential nonnegativity of A with 
respect to K. 

The hit and hold problem considers a positively invariant simplicial cone K in R’, 
and determines the set of initial values X,(K) in R" which eventually reach K (and 
hence remain there due to the positive invariance of K) under x’ = Ax. Neumann and 
Stern (1985) show that X,(K) is itselfa cone, and they characterize it for diagonalizable 
systems A. Berman, Neumann, and Stern (1986) extend these results to 
nondiagonalizable systems A. These results require the spectrum of A to be real, and a 
result from Berman, Neumann, and Stern (1991) requires the matrix A to be essentially 
nonnegative as a condition of the theorem. 

In the case where A is not essentially nonnegative, the hit and hold problem has 
not been solved. However, we use our techniques to find a nonnegative cone K; with 
respect to which A is essentially nonnegative, and the theorem of Berman ez. a/., can be 
applied to a smaller cone K;. If X,(K,) > K, then this is a solution to the hit and hold 
problem for K. 

These results have applications to economics (see, for example, Sierksma, 1979), 


engineering, and biology, where frequently state variables are nonnegative and the 
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essential nonnegativity of the differential system is often required. If a problem is defined 


on the nonnegative orthant and the operator is not essentially nonnegative, our solutions 
may be used to find a sub-cone of R* for which the problem has the required properties, 


and to which the results from control theory may be applied. 

Bitsoris (1988) studied the hit and hold problem for discrete-time dynamical 
systems X;+; = Ax;, and showed that in addition to invariant cones, invariant polyhedral 
sets also arise, and these results could be applied to the unsolved problem of determining 
the spectral properties of systems possessing positively invariant polyhedral cones. For 
the continuous-time case, Castelan and Hennet (1993) again related the existence of 
positively invariant polyhedral sets to the problem of finding an essentially nonnegative 
matrix. 

Since we have a technique for finding essentially nonnegative matrices with respect 
to smaller cones (or larger cones if we consider the transpose of the original system) then 
depending on the requirements of the application we also have a useful technique for 
finding positively invariant polyhedral sets in dynamical systems. 

The requirement that a change of basis be nonnegative is a natural one for many of 
these applications, although the techniques we developed can be extended to cones 
outside of the nonnegative orthant (where for the real spectrum, the problem of seeking an 
essentially nonnegative matrix becomes trivial, and for the general spectrum is simply the 
nonnegative inverse eigenvalue problem for the shifted spectrum). Furthermore, our 
results for complex spectra may be used to extend current results which exist only for 


matrices with real spectra. 
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Stern (1982) shows that when the general cone K = R°, the problem of positively 
invariant cones has its solution in the theory of M-matrices. An /-matrix (after 
Minkowsk1i) is a square matrix A=rl—-B where B20 and r2p(B) (or r> p(B) for 
nonsingular M-matrices). This is a subset of the class of essentially nonpositive matrices 
with nonnegative diagonals. M-matrices occur frequently in problems from the biological 
physical, and social sciences (see, for example, Plemmons, 1977) as well as from the 
mathematical sciences, and Varga (1976) points out that there is an infrequent exchange of 
results between these disciplines. 

Stern (1981) shows that while a matrix which is cone invariant under the 
differential system x’ = Ax, x(0) € K is more general than the requirements for A to be 
an M-matrix. However, the main properties which hold for M-matrices also hold for A, 
which he calls a generalized M-matrix. Numerous classifications of M-matrices have been 
shown (see, for example, Rothblum, 1979; or Stern and Tsatsumeros, 1987) and most of 
these are related to the cone invariance of A or the exponential nonnegativity of —A. 

In terms of other applications of M-matrices, Araki (1975) shows a direct 
application to the method of vector Lyapunov functions (similar to what we presented in 
Chapter II), and Kielson and Styan (1973) relate them to the theory of Markov chains 
which they use to establish results concerning the nonnegative inverse eigenvalue problem 
for complex spectra. 

The negative of an M-matrix is an essentially nonnegative matrix, but the problem 
we present has a less obvious but more important connection to the theory of M-matrices, 


although our solution technique does not provide an immediate solution. The inverse M- 
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matrix problem (see Willoughby, 1977) asks when the inverse of a positive matrix is an M- 
matrix, and it is also mentioned that nonnegative matrices which are not strictly positive 
may also have M-matrices as their inverses, but the zeros may only be those associated 
with reducibility. Since the spectra of M-matrices have certain properties, not every 
positive matrix solves the inverse M-matrix problem, but as our problem allows a shift in 
the spectrum, an irreducible solution to our problem may be a way to obtain an M-matrix 
(or its inverse) which might have use in the numerous applications of M-matrices. 

A more general class of matrices than M-matrices is those with nonpositive off- 
diagonal entries (see Berman, Varga, and Ward, 1978), or essentially nonpositive 
matrices. This class has an important subclass where all of the principal minors are 
positive (see Fiedler and Ptak, 1962) which is also a more general class than M-matrices. 
Again, the relationship of these classes of matrices to essentially nonnegative matrices is 
evident, particularly if we are allowed to shift the spectrum. While the techniques we 
present do not easily adjust to find a matrix with nonpositive off-diagonal elements, it may 
be possible to develop a similar technique for essentially nonpositive matrices. 

Further classes of matrices studied in relation to the problem of when a matrix is 
positive on a cone include the cross-positive matrices, which are those matrices A for 
which x € K and y eK are such that <x, y) = 0, then (x, Ay) >0 fora proper, 
simplicial cone K. These matrices are discussed in Schneider and Vidyasagar (1970) and 
am (1975), and the former draws the relationship to essentially nonnegative matrices. 
This relationship is through the exponential of a matrix, where a matrix A is cross- 


positive on K if and only if e* is positive on K. This is true if and only if A is 


17) 


essentially nonnegative on K. While the paper shows the case of K being the 
nonnegative orthant, the change of bases we constructed can make the result more 
general, either m terms of smaller or larger cones. The class of cross-positive matrices 
includes the copositive matrices, which are symmetric and for which K=R® (see 
Haynsworth and Hoffman, 1969). These matrices are also essentially nonnegative, and in 
addition to the applications previously mentioned, they are important in the field of 
mathematical programming. 

The question of when a linear differential operator is quasimonotone with respect 
to a nonnegative cone is related to the question of when a lmear transformation has an 
invariant cone. Birkoff (1967) studied this problem and used the results to give a 
constructive proof of Perron’s theorems. Following Birkoff, Vandergraft (1968) 
characterized the spectral properties of matrices with invariant cones. These results were 
extended by Barker and Tumer (1973) and were related to the nonnegative verse 
eigenvalue problem. These ideas have been used in algebraic Perron-Frobenms theory 
(see Barker and Schneider, 1975) and have been studied geometrically to extend the ideas 
to non-simplicial cones (see Barker, 1973). An algebraic extension of Lowey and 
Schneider (1975) further characterizes matrices with invariant cones. These techniques 
were applied by Vandergraft (1972) to establish the convergence of some splitting 
methods for solving systems of linear equations. For example, both the Jacobi and Gauss- 
Seidel methods converge for essentially nonpositive matrices with positive inverses (see 
Varga, 1962), although to apply our techniques to this problem may be more challenging 


numerically then solving the original problem Ax = b. 
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The question of when a matrix is quasimonotone nondecreasing (or nonincreasing) 
with respect to a cone is not new, although the requirement that the cone be nonnegative 
does not appeared frequently in recent applications. In applications where the state 
variables are required to be nonnegative (for example, in population dynamics), and when 
the original problem is defined on the nonnegative orthant, this requirement is important. 
These applications, in addition to the one of cone-valued vector Lyapunov functions, 


demonstrate the usefulness of the results of this dissertation. 


a 





VH. CONCLUSIONS AND FUTURE RESEARCH 


This dissertation asks, and answers, the question of when a matnx 
A € R™ is quasimonotone nondecreasing with respect to a nonnegative cone. Chapter 
Il motivates this question with the theory of cone-valued vector Lyapunov functions from 
stability theory in dynamical systems. Chapter VI presents other applications and related 
problems where the results can be used. We summarize these results as follows. 

A matrix A € R"™ is quasimonotone nondecreasing with respect to a nonnegative 
cone if the following conditions hold. 
(1) (Heikkila, 1980). It is necessary that the matrix A have a real first eigenvalue i, 
such that A, = Re(A;) for 1=2,...,n. Furthermore, it is necessary that associated with 
A, be a first eigenvector x; 2 0. 
(2) When o(A), the spectrum of A, is real, it is sufficient that x, > 0. 


A, Ai 


(3) When ofA) CR, x; € (R®?)° © 04 and A= 
O° AS 


is reduced with 


p(Ai1) > p(Az2), it is necessary and sufficient that A2. €¢ R™* be quasimonotone 
nondecreasing with respect to a nonnegative cone. 
(4) When o(A) CR, x; € (R® )° © 04, and A is irreducible (or not reduced) it is 


a 


necessary that for some permutation x, € (R* )° ®0° of x), and for the transformation 


_|() 0 ; al fA, a j ; 
B=]- with A; = B AB = , that the matrix A 22 be quasimonotone 
Xo ys obNes 


nondecreasing with respect to a nonnegative cone. In this case it is sufficient that A; bea 


simple eigenvalue of A. 
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A, A» 


(5) When o(A) CR and a-| 
0 A, 


| satisfies the necessary conditions in (3) 


above, A = p(Aii) = p(Az2), and Aj, and A22 are i upper-triangular quasimonotone 
nondecreasing form with Aj+1,+1 = A, then it is sufficient that Aj,+;20 for 1<i<p. 

(6) For o(A) ZR, it is necessary that x, have positive components corresponding to A, 
and at least all of the eigenvalues with nonnegative imaginary parts. Furthermore, 
considering just the irreducible quasimonotone forms (for which x > 0) it is necessary 
and sufficient that the (shifted) spectrum of A solve the wreducible nonnegative inverse 
eigenvalue problem. 

There are still some unresolved issues for the construction of the cone K(B). The 
case discussed in (5) above (Proposition IV.B.3) is not fully resolved. However, in order 
to apply the construction in (4) above (the proof of Theorem IV.B.2) to such matrices, the 
sufficient condition in (5) is also necessary. 

In Lemma [V.A.1, constructing the matrix B’ with nonnegative inverse B, using 
equation IV.A. 1a requires selecting a basis vector e,. Aside from trying all n 
possibilities, we do not have an efficient method for selecting e. However, this is only 
needed when constructing the cone. The application of vector Lyapunov functionsdoes 
not require the cone’s construction, since it is sufficient to know such a cone exists to 
draw a stability conclusion. 

Furthermore, in Section IV.C we replace some e, k € {1,..., p} with x; when 
constructing the first cone. Again, we may have to try all p possibilities, although we did 


present a method of reducing the search. This is more important to applications that only 
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require the existence of a cone, as this constructive step is required before we can 
determine the quasimonotonicity of A2, and hence, A. 

A problem with the above technique, which is certainly a significant problem in 
applications that require constructing the cone, is that we must compute the generalized 
eigenvectors of a matrix. In case (4) above, if existence is all we have to determine, then 
it may be sufficient to compute the first two eigenvectors of A. However, to find the 
cone we need to compute all of the generalized left eigenvectors of A. While 
theoretically this construction always exists, and while our results give strong existence 
conclusions, it may be numerically intractible to actually construct the required cone. 

Concerning the related problems discussed in Chapter VI, for example the 
applications of M-matrices, in order to use the techniques of this paper to transform a 
given matrix into one of a particular form, then the cone K(B) must certainly be 
constructed. We do not yet know ofa simpler construction than the one we presented. 

When the spectrum is complex we have the further problem of determining which 
part of the real spectrum is required to solve the irreducible nonnegative inverse 
eigenvalue problem. Certainly we can use (for an irreducible, quasimonotone matrix) no 
more eigenvalues than there are nonzero components of x, but we could have a further 
grouping of the complex eigenvalues with different collections of real eigenvalues when 
the final quasimonotone form of the matrix is reducible. We do not know a priori which 
collections may work, aside from trying all of them. As in all of the above problems, since 


the method of vector Lyapunov functions is a proven method in examining large scale 
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dynamical systems (see Siljak, 1978; or Michel and Miller, 1977) we cannot expect to 
have simple numerical results when we employ these techniques. 

Finally, as long as the nonnegative inverse eigenvalue problem is unsolved, then 
this problem is unsolved for the complex spectrum as well, for we have shown that under 
the condition of a positive first eigenvector, the problems are equivalent. 

These questions present further directions for research, if not in terms of the actual 
solution to the problem, then at least in terms of the computability and efficiency of the 
method of solution, particularly when the application requires constructing the cone. We 
have provided analysts and researchers a basic theory and some useful techniques for 


approaching these problems. 
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